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2 R. M. McDermid et al. 

ABSTRACT 

We present the stellar population content of early-type galaxies from the ATLAS^'^ 
survey. Using spectra integrated within apertures covering up to one effective radius, 
we apply two methods: one based on measuring line-strength indices and applying 
single stellar population (SSP) models to derive SSP-equivalent values of stellar age, 
metallicity, and alpha enhancement; and one based on spectral htting to derive non- 
parametric star-formation histories, mass-weighted average values of age, metallicity, 
and half-mass formation timescales. Using homogeneously derived effective radii and 
dynamically-determined galaxy masses, we present the distribution of stellar popula¬ 
tion parameters on the Mass Plane (Mjamj CTe, showing that at fixed mass, 

compact early-type galaxies are on average older, more metal-rich, and more alpha- 
enhanced than their larger counterparts. 

From non-parametric star-format ion histories, we find that the duration of star 
formation is systematically more extended in lower mass objects. Assuming that our 
sample represents most of the stellar content of today’s local Universe, approximately 
50% of all stars formed within the first 2 Gyr following the big bang. Most of these stars 
reside today in the most massive galaxies (> which themselves formed 90% 

of their stars by z ~2. The lower-mass objects, in contrast, have formed barely half their 
stars in this time interval. Stellar population properties are independent of environment 
over two orders of magnitude in local density, varying only with galaxy mass. In the 
highest-density regions of our volume (dominated by the Virgo cluster), galaxies are 
older, alpha-enhanced and have shorter star-format ion histories with respect to lower 
density regions. 

Key words: galaxies: abundances galaxies: elliptical and lenticular, cD galaxies: 
evolution galaxies: stellar content 


1 INTRODUCTION 

Our understanding of the stellar populations in early-type 
galaxies (hereafter ETGs) has made significant progress in 
recent years, through a combination of advancements in both 
observations and stellar population models. These systems 
are often cited as ideal laboratories for uncovering the fossil 
record of galaxy formation processes thanks to the apparent 
lack of dust or star formation, which obfuscates the stellar 
continuum in objects of later type. However, the modern 
observational picture is not quite so straightforward. ETGs 
are now known to host some degree of young stellar popula¬ 
tions (e.g. Trager et al. 2000; Kaviraj et al. 2007; Kuntschner 
et al. 2010). More surprising is the recent finding that siz¬ 
able reservoirs of molecular gas can be found residing inside 
ETGs (Young et al. 2008), coincident with evidence of PAH 
features in the infrared domain (Shapiro et al. 2010), both of 
which strongly indicate ongoing star-formation within these 
galaxies, albeit at a rather low level (Davis et al. 2014). 

In addition to a large range in ages, ETGs are also 
found to show strong correlations between their population 
parameters and mass-sensitive quantities such as luminos¬ 
ity and velocity dispersion (e.g. Jprgensen 1999; Trager et 
al. 2000; Thomas et al. 2005). In general, these studies find 
that high-mass galaxies appear older, metal rich and with en¬ 
hanced abundance ratios with respect to lower mass objects 
- a picture emphatically confirmed with very large samples 
of galaxies from modern surveys (Bernardi et al. 2006; Gal- 
lazzi et al. 2006; Graves et al. 2009a; Thomas et al. 2010). In 
addition to this general behaviour, specific co-dependencies 
between population parameters after controlling for the gen¬ 
eral mass-dependency have also been proposed (Trager et 


al. 2000; Smith et al. 2008; Graves & Faber 2010). Under¬ 
standing what drives the large intrinsic scatter in the stellar 
population scaling relations is key to understanding the evo¬ 
lutionary processes involved. 

Despite a deepening clarity in our observational view of 
ETGs, many open questions remain. For example, the role of 
young stars in these galaxies is still not clear. Although often 
referred to as ‘rejuvenated’ galaxies, it is currently unclear 
if these objects have actually rekindled their star formation 
in recent times, or contrarily, are in the process of recently 
shutting down star formation for the first time. The role of 
environment is also somewhat unclear. Previously consid¬ 
ered as important for ETGs in general (Thomas et al. 2005; 
Schawinski et al. 2007), the interplay of mass and environ¬ 
ment may also play a role (Thomas et al. 2010). 

Here we present the stellar population properties of the 
ATLAS^^ sample (Cappellari et al. 2011a) using two ap¬ 
proaches: one which uses line strengths to derive single stel¬ 
lar population- (SSP-) equivalent parameters to describe the 
age, metallicity and alpha-enhancement; and one which uses 
spectral fitting to obtain a non-parametric star formation 
history that can be used to derive mass-weighted average 
age and metallicity, and estimate the duration of star for¬ 
mation. We compare both approaches, which yield comple¬ 
mentary information, and check for self-consistency between 
them. We then examine the trends of stellar populations with 
various other properties. The ATLAS®^ project comprises 
a well-characterized sample with various multi-wavelength 
data in addition to the SAURON integral field spectroscopy, 
making it a powerful data set to explore correlations between 
multiple properties. 
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This paper proceeds as follows: Section 2 describes 
the observations and analysis methods; Sections 3 and 4 
present our population parameters on the mass-size plane, 
and against velocity dispersion and mass, respectively; Sec¬ 
tion 5 presents trends from our empirical star formation his¬ 
tories, and Section 6 concludes. 


2 OBSERVATIONS, LINE-STRENGTHS AND 

ANALYSIS 

Here we describe the optical spectroscopy used in our analy¬ 
sis. The same spectral data cubes were used for deriving the 
line strengths and for the spectral fitting, so calibration and 
removal of emission applies to both techniques. 

2.1 SAURON observations and basic calibration 

The SAURON spectrograph and basic data reduction is de¬ 
scribed in detail in Bacon et al. (2001). Further details of 
the ATLAS^® SAURON observations are given in Cappellari et 
al. (2011a). The details of our calibrations and line-strength 
measurements follow closely the techniques developed and 
described by Kuntschner et al. (2006). Most of the 260 
ATLAS^^ SAURON spectra were obtained after the spectro¬ 
graph was upgraded with a volume-phase holographic (VPH) 
grating in 2004, except for 64 objects observed as part of 
the SAURON Survey (de Zeeuw et al. 2002) and other related 
projects. The pre- and post-VPH data are reduced and anal¬ 
ysed self-consistently to respect the slightly different disper¬ 
sion characteristics. The wavelength coverage of the pre- and 
post-VPH data are very similar, with only a slight difference 
in resolution. The adopted spectral resolution for pre-VPH 
data is 110 kms“^ compared to 99 kms“^ for VPH data. 
Observations of the almost featureless white dwarf EG131 
were used to verify that the detailed transmission charac¬ 
teristics, including the medium-scale variations described in 
Kuntschner et al. (2006), were present in the pre- and post- 
VPH data, agreeing to better than 1%. Corrections for air- 
mass were applied using the average La Palma extinction 
curve (King 1985). 

A total of 67 different standard stars were observed 
during both the SAURON and ATLAS^*^ campaigns in or¬ 
der to calibrate our line-strengths onto the standard Lick 
system (Worthey et al. 1994), as well as to monitor sta¬ 
bility and reproducibility in our measurements (160 sepa¬ 
rate observations were made). Offsets to the Lick system 
are derived in Appendix A, along with a comparison of pre- 
and post-VPH values, a comparison with the MILES library 
(Sanchez-Blazquez et al. 2006), and discussion of the intrin¬ 
sic uncertainties in our measurements. Appendix B presents 
a comparison of our line-strength measurements on the Lick 
system with Gonzalez (1993) and Trager et al. (1998) for 
overlapping galaxies, showing good agreement given the un¬ 
certainties of matching apertures and differing treatment of 
emission. 

2.2 Emission correction 

The SAURON spectral range includes several nebular emission 
lines, all of which lie close to or within the bandpasses of our 
line-strength indices. Namely: H/3, present in both emission 


and absorption; [0III]AA4959, 5007, which strongly affects 
the Fe5015 index; and [Nl]AA5198, 5200, which falls within 
the Mgb index. As in Kuntschner et al. (2006), we employ 
the technique developed by Sarzi et al. (2006), which in¬ 
volves simultaneously finding the best linear combination of 
absorption- and emission-line templates whilst fitting non- 
linearly for the emission line kinematics. The kinematics of 
the various lines are tied to that of [OIII], being the more 
robustly measured line. The resulting emission-line fits are 
subtracted from the galaxy spectrum if the AoN (ratio of 
the amplitude of the fitted line to the local noise level in the 
spectrum, which includes ht residuals) is significant. Follow¬ 
ing Sarzi et al. (2006), the AoN thresholds used are AoN> 4 
for [OIII], and then if this is detected, AoN> 3 for H/? and 
AoN> 4 for [Nl]. 

The nebular emission lines are fitted and subtracted on 
the spatially-resolved data-cubes (after spatially binning to 
a threshold signal-to-noise ratio of 40 per pixel with the 
Voronoi binning software of Cappellari & Copin 2003) be¬ 
fore generating the aperture spectra. This largely avoids 
strongly non-Gaussian emission-line profiles (due to the ro¬ 
tation and/or distribution of the gas), which can leave sig¬ 
nificant residuals. 

To estimate how strongly the various indices are affected 
by residual under- or over-subtracted emission, we use the 
spectral template fits to provide a good representation of the 
spectral lines in the absence of any emission. We then com¬ 
pare the observed indices with those measured on the tem¬ 
plate fits. The differences between the observed and template 
indices are typically similar to or smaller than the statisti¬ 
cal errors on the line measurements. However, a few galaxies 
show significantly larger differences, due to imperfect emis¬ 
sion correction in the presence of very strong emission lines 
with complex, non-Gaussian line profiles. Rather than use 
the index measured on the template, which is known to be 
biased toward solar abundance ratios, we adopt the differ¬ 
ence between the observed and best-fit template indices as an 
estimate of the systematic errors. We furthermore take the 
statistical uncertainties from repeat observations of standard 
stars as a lower limit on the errors (see Table A2 given in 
Appendix A). 

2.3 Aperture line strength measurements 

Figure 1 of Emsellem et al. (2011) shows the distribution of 
radii covered by our integral-field spectroscopic data, defined 
as the maximum circular radius, Rmax, that is at least 85% 
filled with spectra. Note that this fill factor includes lost 
spaxels removed due to e.g. foreground stars. The median 
maximal radius is ~1.2 Re, with only 10 galaxies having 
coverage of less than Rel2. In terms of aperture coverage, 
more than half the sample has a filling factor of 85% or 
more within a circle of Re, and all but 11 galaxies have this 
coverage at Rel‘1- 

We measure line strengths in three different circular 
aperture sizes: Re!^, Re/2 and Re- The smallest of these is 
comparable to aperture sizes commonly used in conventional 
long-slit investigations (e.g. Trager et al. 2000; Thomas et al. 
2005). The Re/2 aperture is our largest common aperture 
that requires minimal aperture corrections in a small num¬ 
ber of objects. The Re aperture is the most representative 
‘global’ value that we can derive based on our data. 
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Figure 1. Absorption line indices measured within Re/S against the effective velocity dispersion, cje (measured within one Re, taken 
from Cappellari et al. 2013a). Red symbols and text indicate problematic galaxies as described in the main text. The dotted magenta 
line in the Fe5015 plot indicates the separation of ‘low-Fe’ objects discussed in the text, given here as log(Fe5015) = 0.53. The solid line 
in the bottom panel indicates the biweight mean value, 1.84, with dashed lines showing the ±5 times the standard deviation around this 
value. Closed and open symbols in the upper panels correspond to galaxies with and without Fe5270 measurements respectively. 
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Figure 2. Index-index diagram for our sample using the indices 
H/3 versus [MgFe50]', overplotted with model predictions from 
Schiavon (2007) for single stellar populations with constant ages, 
metallicities and solar abundance ratios. From top to bottom, each 
panel uses indices measured within different circular apertures, 
with radii of Re/S, i?e/2 and Re as indicated. An error bar in 
each panel indicates the statistical errors. 


To generate the aperture measurements from the data- 
cubes, we use an analogy of ‘growth curves’, co-adding spec¬ 
tra within increasing circular aperture sizes and measuring 
the indices at each step. We assign a characteristic radius to 
each aperture spectrum that represents the radius of a circle 
having the same area as that of the included SAURQN data, 
computed from the actual bins within the aperture. Thus, 
where the aperture is well-filled by the SAURON data, this ra¬ 
dius is similar to the radius of the desired circular aperture. 
But as the aperture becomes more poorly filled near the edge 
of the field, the characteristic radius will tend to a hxed value 
(the area of the full observed field) regardless of the target 
aperture size. By using a threshold filling factor of 85%, the 
difference between the target aperture and delivered area is 
small. 

For the galaxies where the maximum characteristic ra¬ 
dius is smaller than the target aperture radius, an aperture 
correction is applied to the measured indices. Due to the 


presence of age gradients that systematically vary in strength 
with galaxy age, we derive age-dependent aperture correc¬ 
tions. Appendix C gives details of this correction and of the 
coefficients to apply such corrections. The aperture correc¬ 
tions for our sample are typically 2% or less, with a maxi¬ 
mum of 8% for our largest correction, and our overall results 
do not depend on these corrections. 

The line indices themselves are measured on the 
Lick/IDS system (Worthey et al. 1994) following the stan¬ 
dard treatment of broadening the spectra as per the Lick 
broadening function. We account for the dilution of the 
line strength due to velocity broadening by measuring the 
equivalent indices on the broadened and unbroadened op¬ 
timal template fit, thus accounting for line-strength and 
LOSVD dependencies of the velocity dispersion corrections 
(Kuntschner et al. 2006; McDermid et al. 2006). We measure 
the following indices in all galaxies and aperture sizes: H/3, 
Fe5015 and Mg&. The SAURQN wavelength range has a red 
cut-off wavelength that changes as a function of position in 
the field, causing Kuntschner et al. (2006) to dehne a new 
index, Fe5270s, which linearly traces the Fe5270 Lick index. 
We use this modihed index to estimate the standard Lick 
Fe5270 index for a subset of objects and apertures where the 
SAURON wavelength coverage allows.^ The line strengths mea¬ 
sured within the three aperture sizes considered are given in 
Tables 1-3, which are available in full via the online journal, 
and also via our project webite http://purl.com/atlas3d. 

Figure 1 presents our line indices measured within one 
of the three apertures. Re/8, plotted as a function of the 
velocity dispersion measured within one Re, log((Je) taken 
from Cappellari et al. (2013a). Filled symbols denote the 
galaxies where Fe5270 could be measured within the given 
aperture size, with the remainder plotted as open symbols. 

As found by numerous other authors (e.g. Trager et al. 
2000; Thomas et al. 2005), there are rather tight, positive 
correlations between the metal-sensitive lines and velocity 
dispersion, and the age-sensitive H/3 index shows a negative 
trend with Ue. All indices show an increased scatter at lower 
CTe values, implying a greater variety of apparent stellar pop¬ 
ulation properties at lower masses. 

The Fe5015 index lies on top of the strong [OIII] emis¬ 
sion line at 5007A. For this index, the apparent scatter is 
driven by a number of galaxies sitting below the main rela¬ 
tion. To test whether this is the result of infill from poorly- 
subtracted [OIII] emission or genuinely weaker iron, we use 
the fact that some galaxies have a measurable Fe5270 index, 
which is free from nebular emission contamination. Figure 
1 shows several galaxies where Fe5270 also lies below the 
general relation, but the ratio of the two iron indices (bot¬ 
tom panel of Figure 1) is consistent with the sample average, 
demonstrating that the weak Fe5015 measurements in these 
cases are not due to poor emission-line correction. We discuss 
these ‘low-Fe’ objects further in Section 2.6. 

We highlight a few objects in Figure 1 that need special 
treatment. Firstly, two objects had signihcantly lower signal- 
to-noise ratio than the rest of the sample, namely NGC4268 
and PGC 170172. These two objects are not considered fur¬ 
ther in our analysis. In addition, the following objects show 


^ We derive Fe5270 using the empirical relation Fe5270 = 1.28 X 
Fe5270s + 0.03 from Kuntschner et al. (2006) 
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large errors in their Fe5015 values due to imperfect subtrac¬ 
tion of the very strong and non-Gaussian [OIII] emission 
line, as indicated by their spuriously low Fe5015/Fe5270 ra¬ 
tio in Figure 1: NGC 1222 (starburst galaxy, Balzano 1983), 
NGC4684, NGC 5273, and NGC 3998. These latter two ob¬ 
jects are known LINERs with a broad nuclear emission line 
component (Veron-Cetty & Veron 2010) that only affects the 
small Re/S aperture size, and so we include them when con¬ 
sidering the larger aperture measurements. NGC 1222 is ex¬ 
cluded for all aperture sizes. PGC 029321 has no Fe5270 mea¬ 
surement, but shows large discrepancies between the Fe5015 
index measured on the observed and fitted spectrum due to 
emission residuals, and so is also excluded in our analysis. 
For all the problematic cases mentioned here, we still pro¬ 
vide measured values in Tables 1-3 where possible, but flag 
them in the ‘quality’ column. 

2.4 Deriving SSP parameters 

As in McDermid et al. (2006), we make use of a customized 
version of the Schiavon (2007) models where [Ti/Fe] = 0.5 x 
[Mg/Fe]. These models were found to give better agreement 
between abundance ratios computed with Fe5015 and Fe5270 
(see Kuntschner et al. 2010, for details), and this relative 
scaling of [Ti/Fe] and [Mg/Fe] is also approximately consis¬ 
tent with the measured abundances for early-type galaxies 
from Johansson et al. (2012). The models make predictions 
of the Lick indices for a grid of various ages, metallicities and 
alpha-element enhancements. Each position in this grid rep¬ 
resents the prediction for a single population of stars sharing 
these parameters. 

We show in Figure 2 the distribution of our full sam¬ 
ple in the indices H/l versus [MgFe50]'^ for the three aper¬ 
ture sizes, overplotted with these model predictions using 
[a/Fe[= 0. The effects of abundance ratio are largely miti¬ 
gated with this choice of indices (Thomas et al. 2003). The 
models have been adapted for zero-point differences between 
the Lick system and the Schiavon models using Table 1 of 
Schiavon (2007). 

From this figure, it can be seen how the changing dis¬ 
tribution of H/3 for the different apertures implies that the 
young stars in general are centrally concentrated. It can also 
be seen that the inferred metallicity becomes lower for larger 
apertures, due to the inclusion of the metal-poor outer re¬ 
gions. We also note that several objects have implied ages 
that are older than the canonical age of the Universe from 
Planck (13.798 ±0.037 Gyr; Planck Gollaboration, Ade et al. 
2013). In Appendix D we show that these points are consis¬ 
tent with objects close to this upper age limit considering 
our observational uncertainties. 

To derive the best-fitting SSP parameters for our data, 
we use the three indices that are measured across the full 
field for all galaxies, namely H/3, Fe5015 and Mgb. We em¬ 
ploy the straightforward yet robust technique of finding the 
SSP model that best predicts our observed line indices simul¬ 
taneously. As per McDermid et al. (2006) and Kuntschner et 
al. (2010), we do this by means of chi-squared fitting: effec¬ 
tively computing the ‘distance’ from our measured indices 
to all the predicted values of those indices, and finding the 

2 [MgFeSO]' = 0.69xMgt,+Fe50l5 (Kuntschner et al. 2010) 


model with the minimum total distance. To reduce the effects 
of grid discretization, we oversample the original models us¬ 
ing linear interpolation. Measurement errors on the indices 
are included as weighting factors in the sum. This particular 
technique has been used by several authors (Proctor et al. 
2004; Thomas et al. 2010), although there are other methods 
for inverting the model grids to derive parameter values (e.g. 
Graves & Schiavon 2008). 

We estimate the parameter errors by taking the range 
of points at which the difference in chi-square from the min¬ 
imum is equal to unity, giving an estimate of the projected 
dispersion on each parameter. This method is computation¬ 
ally efficient, and we verified that it agrees well with a Monte- 
Carlo approach. All SSP values are given in Tables 1-3, 
which are available in full via the online journal, and via 
our project website: http://purl.com/atlas3d.. 

2.5 Mass-weighted parameters from spectral 
fitting 

We use the publicly-available^ penalized pixel fitting code 
(pPXF, Cappellari & Emsellem 2004) to fit a linear com¬ 
bination of SSP model spectra from the MIUSCAT model 
library (Vazdekis et al. 2012), spanning a regular grid of 
log(age) (0.1 - 14 Gyr) and metallicity {[Z/H] = —1.71 — 
0.22, equivalent to Z = 0.0004 — 0.03), giving 264 templates 
in total. Due to the degeneracy of inferring the star forma¬ 
tion history from the integrated galaxy light, we make use of 
linear regularization (Press et al. 1992) to impose a smooth¬ 
ness constraint on the solution, using the ‘REGUL’ option of 
pPXF. The solution is defined by the weights applied to each 
model in a regularly-sampled grid of age and metallicity, and 
regularisation imposes a constraint in this parameter-space 
such that models with neighbouring age and metallicity must 
have weights that vary smoothly. The relative importance of 
this smoothness constraint compared to the actual fit quality 
is controlled via a simple ‘regularisation parameter’, which 
essentially controls how smooth the solution is. The SSP 
model spectra are provided to an initial birth cloud mass 
of 1 solar mass, and so the relative ‘zero-age’ mass-to-light 
ratios are included in the fitted weights. The distribution of 
weights from pPXF therefore defines the relative contribu¬ 
tions in (zero-age) mass of the different populations: effec¬ 
tively the star formation history (SFH) for that galaxy. 

For each galaxy, we fit the spectrum integrated within 
one effective radius as per the largest aperture considered 
for the SSP analysis. We derive solutions for a range of reg¬ 
ularization parameter values, and select the value that gives 
a resulting best-fit solution with Ax^ ~ y/^N, where N is 
the number of pixels in the spectrum being fitted, and Ax^ 
is computed as the difference between the x^ of tho current 
solution with that of the non-regularized case. This criterion 
is suggested by Press et al. (1992) as a robust indication of 
when the fit has become unacceptable. In this way, we are 
deriving the smoothest star formation history that still re¬ 
produces the SAURON spectrum. The actual star formation 
histories of galaxies are likely not smooth. Star formation is 
inherently stochastic when considered with sufficient tempo¬ 
ral and spatial resolution, and specific events such as mergers 

^ http://purl.org/cappellari/idl 
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Table 1. Measured Lick index measurements and SSP parameters measured within Re/S. 


Name 

(1) 

H/3 [A] 

Re/S 

(2) 

Fe5015 [A] 

Re/S 

(3) 

Mg 6 [A] 

Re/S 

(4) 

Fe5270 [A] 

Re/S 

(5) 

Agessp [Gyr] 
Re/S 
(6) 

[Z/H]ssp 

Re/S 

(7) 

[o/FeJssp 

Re/S 

(8) 

Quality 

(9) 

IC0560 

2.95 ± 0.10 

4.17 ± 0.14 

2.29 ± 0.12 

2.20 ± 0.11 

1.85 ± 0.21 

-O.lli 0.05 

0.09 ± 0.05 

1 

IC0598 

2.88 ± 0.10 

4.64 ± 0.14 

2.46 ± 0.12 

2.31 ± 0.12 

1.70 ± 0.22 

0.04 ± 0.05 

0.05 ± 0.04 

1 

IC0676 

3.51 ± 0.15 

2.64 ± 0.14 

1.75 ± 0.23 

1.51 ± 0.11 

2.02 ± 0.23 

-0.64 ± 0.08 

0.37 ± 0.11 

2 

IC0719 

2.16 ± 0.10 

5.14 ± 0.14 

3.26 ± 0.12 

2.75 ± 0.11 

3.53 ± 0.68 

0.03 ± 0.05 

0.11 ± 0.04 

1 

IC0782 

2.27 ± 0.10 

4.98 ± 0.14 

2.88 ± 0.12 

- ± - 

3.10 ± 0.36 

-0.03 ± 0.06 

0.06 ± 0.05 

1 

IC1024 

4.00 ± 0.10 

2.35 ± 0.14 

1.52 ± 0.17 

1.20 ± 0.11 

1.59 ± 0.18 

-0.75 ± 0.06 

0.41 ± 0.09 

2 


Columns: (1) Galaxy name. (2)-(5) Lick index measurements within Re/S, measured on the Lick/IDS system. Fe5270 is measured 
indirectly via the Fe5270s by applying the relation Fe5270 = 1.28 X Fe5270s + 0.03 (Kuntschner et al. 2006). Due to field coverage 
issues, this index cannot be measured for all objects. (6)-(8) Estimates of SSP-equivalent age, metallicity and abundance ratio within a 
circular aperture of radius Re/S using the SSP models of Schiavon (2007) and using the three indices H/3, FeSOlS, and Mgfe. (9) Quality 
flag: 1 = Data are of good quality; 2 = Object with weak FeSOlS index in the Re/S aperture; 3 = Strong emission line residuals; 4 = 
low signal-to-noise ratio. The full version of this table with all 260 galaxies is made available via the online journal, and via our project 

website http://purl.com/atlas3d. 


Table 2. Measured Lick index measurements and SSP parameters measured within i7e/2. 


Name 

(1) 

H/3 [A] 

Re/2 

(2) 

Fe5015 [A] 
Re/2 

(3) 

Mg 5 [A] 

Re/2 

(4) 

Fe5270 [A] 
Re/2 

(5) 

Agessp [Gyr] 
Re/2 
(6) 

[Z/H]ssp 

Re/2 

(7) 

[o/Fe]ssp 

Re/2 

(8) 

Quality 

(9) 

IC0560 

2.53 ± 0.10 

3.86 ± 0.14 

2.32 ± 0.12 

2.15 ± 0.11 

2.67 ± 0.31 

-0.30 ± 0.05 

0.14 ± 0.06 

1 

IC0598 

2.43 ± 0.10 

4.20 ± 0.14 

2.47 ± 0.12 

2.22 ± 0.11 

2.84 ± 0.33 

-0.22 ± 0.05 

0.11 ± 0.05 

1 

IC0676 

3.02 ± 0.10 

3.25 ± 0.14 

1.92 ± 0.12 

1.75 ± 0.11 

2.20 ± 0.25 

-0.48 ± 0.06 

0.22 ± 0.07 

2 

IC0719 

2.43 ± 0.10 

4.16 ± 0.14 

2.69 ± 0.12 

2.46 ± 0.11 

2.84 ± 0.33 

-0.18 ± 0.04 

0.18 ± 0.05 

1 

IC0782 

2.08 ± 0.10 

4.51 ± 0.20 

2.58 ± 0.12 

- ± - 

4.77 ± 0.85 

-0.25 ± 0.05 

0.03 ± 0.06 

1 

IG1024 

3.59 ± 0.15 

2.25 ± 0.14 

1.66 ± 0.20 

1.30 ± 0.11 

2.06 ± 0.24 

-0.78 ± 0.07 

0.50 ± 0.05 

2 


Columns are as per Table 1 but for an i?e/2 aperture size. The full version of this table with all 260 galaxies is made available via the 
online journal, and via our project website http://purl.com/atlas3d. 


Table 3. Measured Lick index measurements and SSP parameters measured within Re- 


Name 

(1) 

H/3 [A] 

Re 

(2) 

Fe5015 [A] 
Re 
(3) 

Mg 5 [A] 

Re 

(4) 

Fe5270 [A] 
Re 
(5) 

Agessp [Gyr] 
Re 
(6) 

[Z/H]ssp 

Re 

(7) 

[n/Fe]sSP 

Re 

(8) 

Quality 

(9) 

IC0560 

2.37 ± 0.10 

3.75 ± 0.16 

2.31 ± 0.12 

2.10 ± 0.10 

3.53 ± 0.58 

-0.38 ± 0.06 

0.19 ± 0.07 

1 

IC0598 

2.30 ± 0.10 

4.03 ± 0.14 

2.41 ± 0.12 

2.16 ± 0.11 

3.38 ± 0.57 

-0.30 ± 0.05 

0.12 ± 0.05 

1 

IC0676 

2.66 ± 0.10 

3.37 ± 0.14 

1.96 ± 0.13 

1.79 ± 0.11 

2.84 ± 0.38 

-0.51 ± 0.06 

0.21 ± 0.07 

2 

IC0719 

2.35 ± 0.10 

3.84 ± 0.14 

2.56 ± 0.12 

2.33 ± 0.10 

3.38 ± 0.60 

-0.29 ± 0.05 

0.22 ± 0.06 

1 

IC0782 

2.03 ± 0.10 

4.47 ± 0.27 

2.50 ± 0.12 

- ± - 

5.31 ± 0.85 

-0.30 ± 0.06 

0.01 ± 0.07 

1 

IC1024 

3.33 ± 0.18 

2.47 ± 0.14 

1.71 ± 0.16 

1.39 ± 0.12 

2.29 ± 0.42 

-0.73 ± 0.06 

0.43 ± 0.07 

2 


Columns are as per Table 1 but for an Re aperture size. The full version of this table with all 260 galaxies is made available via the 

online journal, and via our project website http://purl.com/atlas3d. 


and interactions may trigger or truncate the star formation 
process on short timescales. In fact, we have also shown that 
the colors and H/3 indices of the bluest and most H2-rich 
galaxies can be reproduced with a small amount of current 
or recent star formation activity superposed on an old stel¬ 
lar population (Young et al. 2014). However, inferring such 
events from the integrated light is highly degenerate for typ¬ 
ical observational data, and the smoothness constraint offers 
a robust and objective result that can be used to study sys¬ 
tematic trends in the entire sample. We also note that non¬ 
smooth solutions are in no way excluded by our approach, 


but if a smooth solution can give a fit of similar quality 
(as determined via Ay^), then the smooth solution will be 
adopted. 

The mass-weighted age and metallicity are given as: 


log(AgesFH) = 


'^Wi log(tssp,i) 


and 


[2/H]sfh ^ 




( 1 ) 

( 2 ) 
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Table 4. Mass-weighted stellar population properties measured 
within Re. The full version of this table with all 260 galaxies is 
made available via the online journal, and via our project website 
http://purl.com/atlasSd. 


Name 

(1) 

AgespH [Gyr] 
Re 
(2) 

[Z/H]sfh 

Re 

(3) 

Go [Gyr] 

Re 

(4) 

IC0560 

6.73 ± 1.07 

-0.46 ± 0.06 

6.55 ± 1.07 

IC0598 

5.37 ± 0.99 

-0.40 ± 0.05 

8.60 ± 0.99 

IC0676 

6.00 ± 1.92 

-0.55 ± 0.18 

7.22 ± 1.92 

IC0719 

9.21 ± 0.86 

-0.25 ± 0.08 

4.23 ± 0.86 

IC0782 

6.16 ± 0.83 

-0.19 ± 0.11 

7.09 ± 0.83 

IC1024 

4.39 ± 0.63 

-0.77 ± 0.15 

8.65 ± 0.63 


Columns: (1) Galaxy name. (2) Mass-weighted age. (3) Mass- 
weighted metallicity. (4) Half-mass formation time. These prop¬ 
erties were derived from a regularised spectral fit to a spectrum 
integrated within Re. Errors were computed via Monte-Carlo re¬ 
alisations without regularisation constraints. See text for details. 

respectively, where Wi denotes the weight given by pPXF 
to the i’th template, which has age tssp,i and metallicity 
[Z/H]ssp,i- The grid of model ages is sampled logarithmi¬ 
cally, hence the logarithmic form of equation 1. Errors on 
these quantities were computed via Monte-Carlo simulations 
without regularisation. 

The SSP spectral templates we use have fixed (solar) 
abundance ratios, and so we cannot derive a mass-weighted 
[a/Fe]. We can, however, derive a measure of the star for¬ 
mation timescale directly from the SFH itself. The measure 
we use is the time, tso, taken to form 50% of the current day 
stellar mass within one Re, and is derived from the cumula¬ 
tive function of the SFH (see section 5.1). All mass-weighted 
quantities are given in Table 4. 

2.6 Comparison of SFHs with SSPs and SDSS 

An obvious concern in deriving the SFH from SAURON spec¬ 
tra is the short wavelength range available. To test this, we 
compare the SFHs derived from SAURON with those derived 
from fitting SDSS spectroscopy (DR8, York et al. 2000) of 
the galaxies in common, using the entire SDSS wavelength 
range to constrain the SFH - a more than ten-fold increase 
in the wavelength coverage of the spectrum. For this test, 
we adapt the SAURON aperture to match the 3” diameter 
SDSS fiber, thus minimizing the effects of population gradi¬ 
ents. Figure 3 shows a comparison of the SAURON and SDSS 
mass-weighted average ages and metallicities derived by ap¬ 
plying the above procedure. The mass-weighted parameters 
compare very well, showing no significant systematic biases 
across the range of age and metallicity considered despite 
the dramatically different spectral ranges used. Age shows 
a slight trend to older ages derived by SAURON, but the old¬ 
est case is still consistent with SDSS within the la scatter. 
Metallicity is also very consistent, albeit with a small av¬ 
erage offset, but again within the Icr scatter. Overall, the 
agreement is remarkably good, and suggests that the short 
wavelength range of SAURON does not significantly bias our 
derived SFHs. 

We compare our mass-weighted and SSP-equivalent pa¬ 
rameters in Figure 4, this time using a IRe aperture. The 


A= 0.47 ± 1.17 Gyr 



0 2 4 6 8 10 12 

AgesFH SAURON (Gyr) 

A=-0.05± 0.06 



-0.4 -0.2 0.0 0.2 

[Z/H]sph SAURON 


Figure 3. Comparison of mass-weighted age (top) and metallicity 
(bottom) derived from fitting SAURON and SDSS spectra within the 
same spatial aperture. The mean difference (SAURQN-SDSS) and 
the standard deviation are given in the plot titles. The dotted 
lines indicate the 1;1 identity line, with a thick line showing the 
best linear fit. Errors were computed via Monte-Carlo simulations 
without regularisation. 


mass-weighted ages are systematically older for all but the 
oldest ages, where we are biased by the lower maximum age 
of our spectral template library. Including templates up to 
17 Gyr reduces the ‘bunching’ of points at oldest ages, but for 
clarity of interpreting our SFHs, we use a maximum age of 
14 Gyr in our template library, consistent with current cos¬ 
mological models (Planck Gollaboration, Ade et al. 2013). 
The older mass-weighted ages derived from spectral fitting 
are a direct result of the well-known bias of SSP-equivalent 
ages by small mass fractions of young populations, and our 
findings agree well with the exploration of this issue by Serra 
& Trager (2007). Metallicities show good agreement, also as 
expected following the findings of Serra & Trager (2007), 
who showed that SSP-equivalent metallicities are close to 
the mass-weighted values, being less biased by the presence 
of young stars. We also note that there is no obvious bias 
caused by using models limited to solar abundance ratios. 

Finally, we find a clear trend between the SSP- 
equivalent alpha enhancement and tso. We over-plot the lin¬ 
ear relation given by de La Rosa et al. (2011), who conducted 
a similar analysis on SDSS spectra using the STARLIGHT 
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[^/^]ssp [®/F^®]sSP 


Figure 4. Comparison of SSP-equivalent parameters (x-axis) and mass-weighted parameters derived using spectral fitting (y-axis) within 
the same aperture of one effective radius R^. Magenta points highlight the objects with low Fe as discussed in Section 2.3. Identity lines 
are given for age and metallicity. The right panel shows [a/Fe] versus tso, the time taken to form half the current-day mass. The relation 
iso = —15.3[o/Fe] -|- 5.2 given by eqn. 2 in de La Rosa et al. (2011) is over plotted on our data with a green dashed line. The blue dotted 
line shows the canonical conversion of abundance ratio to star formation timescale from equation (4) of Thomas et al. (2005) based on 
chemical evolution models: [a/Fe] ~ | ~ g log At. The red solid line and text indicate our best-fitting function of this same form, giving 
an empirical relation between [a/Fe] and fso for the ATLAS®*^ sample. 





M(H2)/Mj,„ M(H,)/Mj^„ M(H2)/Mj,„ 

Figure 5. Residuals from the SSP-velocity dispersion trends plotted as a function of molecular gas fraction. Median errors are indicated 
in the bottom left corner of each panel, with H 2 mass fraction errors within the plotting symbol. Arrows indicate 3 (t upper-limits for 
the H 2 mass. Diamonds denote CO detections in the Virgo cluster. Open circles show detections of non-Virgo objects. Filled magenta 
symbols correspond to the ‘low-Fe’ objects discussed in Section 2.3. Their weak Fe5015 indices result in low metallicity and high [o/Fe]. 
We speculate that these objects have had recent gas accretions, giving rise to young, metal poor stars that have not yet been enriched 
with ejecta from Type la supernovae. Only one of these low-Fe objects is found in the Virgo cluster (filled magenta diamond). 


code of Cid Fernandes et al. (2005). We also show the 
canonical relation between abundance ratio and star for¬ 
mation timescale from Thomas et al. (2005), equation (4): 
[a/Fe\ ~ I — I log At, where At is the FWHM of the 
assumed Gaussian SFH. This relation reproduces the non¬ 
linearity of our measured timescale against [a/Fe]. We also 
note that the limited age resolution of the SSP models at the 
oldest ages creates a ‘floor’ in the tso values at the shortest 
timescales. This important region, corresponding to the ear¬ 
liest epochs of massive galaxy formation, is where abundance 
ratios can yield higher fidelity than our SFH. We therefore 
use our empirical half-mass formation time to constrain a 
relation of the same form as that of Thomas et al. (2005) to 


relate [a/Fe] to tso. The resulting expression is indicated by 
the red line in Fig. 4, and is given by: 

[a/Fe] = 0.28 - 0.19 log/tso) • (3) 

The objects plotted with magenta symbols in Figure 
4 (the low-Fe objects) are mostly outliers to the [a/Fe]- 
tso trend, having large values of tso - i-®- long formation 
timescales, contrary to what is expected from their high 
[a/Fe] values. The relatively weak iron indices result in low 
metallicities and subsequently high [a/Fe], but since they 
also have low velocity dispersions, they are strong outliers 
in the scaling relations shown in Figure 8. Figure 5 shows 
the residuals of the SSP parameter - velocity dispersion re¬ 
lations for an Re/8 aperture, plotted against molecular gas 
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mass fraction, with the low-Fe objects marked with filled 
symbols. This shows that, as well as having low metallic- 
ity, the low-Fe objects also have high molecular gas fractions 
and are among the youngest in the sample. This suggests 
their apparently high SSP-equivalent [a/Fe] may indicate 
the presence of stars that have not yet been appreciably 
enriched by Type la supernovae - the main contributor of 
iron to the gas involved in star-formation. Their young age 
causes them to outshine the older, less a-enhanced stars, 
and dominate the SSP abundance measurement. They do 
not show any preferred kinematic misalignment (Davis et al. 
2013), and are not Virgo cluster members except for one, 
NGC4694, which is a low-mass object in the outskirts of 
Virgo. We speculate that these low-Fe galaxies have recently 
accreted a new supply of cold, low-metallicity gas, which is 
currently in the process of forming stars. These twelve low-Fe 
objects are shown as magenta symbols in subsequent plots 
unless stated otherwise, and are indicated in Tables 1-3 via 
the ‘quality’ flag. 


3 POPULATION PARAMETERS ON THE 

MASS PLANE 

Galaxy stellar populations have been found to correlate bet¬ 
ter with stellar surface density (Kauffmann et al. 2003; Franx 
et al. 2008) or with velocity dispersion (Graves & Faber 
2010) than with mass or galaxy size. The same was found to 
hold for the dynamical M/L, which mainly traces age vari¬ 
ations (Gappellari et al. 2006). In Cappellari et al. (2013b) 
we used accurate masses and velocity dispersions to show 
that indeed the velocity dispersion is the best simple tracer 
of multiple galaxy properties. We showed that velocity dis¬ 
persion is such a good predictor of galaxy properties because 
it traces the bulge fraction, which appears to be the main 
driver of the M/L, H/3 and colour, as well as the molecular 
gas fraction in ETGs. Here we focus on trends in the stellar 
population parameters. 

Figure 6 and 7 show the distribution of population pa¬ 
rameters as a function of mass and size (which we hereafter 
refer to as the Mass Plane) for our sample, using the SSP- 
equivalent and mass-weighted parameters respectively. The 
mass Mjam is derived applying Equation 3 of Gappellari et 
al. (2013b), using the mass-to-light ratio (M/L)jam and r- 
band luminosity given in Table 1 of Gappellari et al. (2013a), 
columns 6 and 15 respectively. We use the robust size param¬ 
eter corresponding to the major axis of the half-light 

isophote (Hopkins et al. 2010) as derived in Gappellari et al. 
(2013a), and taken from Table 1 there. 

The points in Figs. 6 and 7 are coloured according to the 
population parameters as indicated. The distributions have 
been adaptively averaged by applying an implementation‘s 
of the LOESS algorithm (Gleveland & Devlin 1988), using a 
local linear approximation, and a regularization factor / = 
0.6. The same approach is taken in Cappellari et al. (2013b) 
and is found to give a good approximation to the underlying 
intrinsic distribution, mimicking the average behaviour that 
would be derived from much larger samples (e.g. Gallazzi et 
al. 2006). We note that, since we are adaptively averaging 

■s available from http://purl.org/cappellari/idl 


the data points, the range of values is reduced in accordance 
with the reduced effects of errors and intrinsic scatter. 

As in Cappellari et al. (2013b), we over-plot lines of con¬ 
stant velocity dispersion, as implied by the virial mass esti¬ 
mator Mjam= 57?™‘*V^/G. We find that the SSP-equivalent 
population parameters follow trends of nearly constant ve¬ 
locity dispersion on the Mass Plane, as indicated by the tilt 
of the iso-colour regions. The mass-weighted age and for¬ 
mation timescale follow even more tightly the lines of con¬ 
stant dispersion, which reflects the lower scatter in the mass- 
weighted properties as a result of lower formal uncertainties, 
and less susceptibility to small fractions of young stars be¬ 
ing present. The mass-weighted metallicity is also closely 
linked to the velocity dispersion, though the iso-colour re¬ 
gions show an appreciable dependence on Mjam also (i.e. 
steeper iso-colour regions). 

These figures show that, for fixed mass, more com¬ 
pact galaxies are on average older, more metal rich, have 
higher abundance ratios and shorter formation timescales 
than larger galaxies. The trends with SSP age and metal¬ 
licity agree with previous results from averaged properties 
of stacked SDSS spectra (van der Wei et al. 2009; Shankar 
et al. 2010). We show here that these trends hold for mea¬ 
surements of individual galaxies. Moreover, the trend with 
abundance ratio, mass-weighted properties, and formation 
timescale measured directly from the star-formation history 
can only be obtained from a spectroscopic analysis, and are 
reported here for the first time. 

Graves et al. (2009b) found that Re does not correlate 
with either population parameters or their residuals from 
the population-(T relations. This can be understood from the 
mass-size plane. At each there is a similar range of 

population parameters, thus showing no particular correla¬ 
tion. Likewise, as a traces the parameters closely, the resid¬ 
uals of the population-(j trends show no significant residual 
correlation with R^^K Indeed, Figures 6 and 7 show that ve¬ 
locity dispersion is an excellent predictor of all three stellar 
population parameters, with the smoothed population pa¬ 
rameters varying essentially parallel to the lines of constant 
(virial-predicted) velocity dispersion. 

We advocate that the Mass Plane is a novel and use¬ 
ful tool for understanding galaxy properties. Gomparing our 
Figures 6 and 7 to Figures 3-8 of Gappellari et al. (2013b) 
reveals the close connections between the stellar population 
parameters and other important quantities, in particular 
colour, molecular gas mass fraction, and bulge fraction, such 
that, at a given mass, smaller galaxies are redder in colour, 
have lower fractions of molecular gas, and have a second 
moment Kms = dominated by a dynamically hot 

central ‘bulge’ component. There are no particular trends 
with intrinsic shape other than the fact that the population 
of round objects above Mjam> 2 x 10^^ M© are typically old, 
metal rich and a—enhanced. Nor is there a close connection 
to the specific angular momentum, Ar, which supports the 
findings of Naab et al. (2014) where Ar is shown to trace 
a galaxy’s evolution, but that there are many evolutionary 
paths to given current-day value. We explore the connec¬ 
tion between the kinematic classes and stellar populations 
in Section 4.5. 
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Figure 6. SSP-equivalent population parameters measured within one effective radius, plotted on the plane of effective radius 
versus dynamical mass, Mjam- Colours indicate the population parameters as indicated in each plot after spatially averaging with the 
LOESS technique described in the text. For a fixed range in Mjam? the smallest galaxies are generally the oldest, most metal rich and 
alpha-enhanced. Dashed lines show lines of constant velocity dispersion: 50, 100, 200, 300, 400 and 500 kms“t from left to right, as 
implied by the virial mass estimator Mjam= 5R^^^cr^/G. The red curve shows the ‘zone of exclusion’ (ZOE) defined in Cappellari et al. 
{2013b). 
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Figure 7. Same as Figure 6, but using the mass-weighted age and metallicity, and tgo - the time taken to form half the stars inside IRe- 


4 POPULATION TRENDS WITH VELOCITY 
DISPERSION AND MASS 

In this section we present a more traditional two-dimensional 
analysis of the stellar population parameters, examining 
their scaling properties with velocity dispersion and mass, 
for SSP-equivalent and mass-weighted properties. This al¬ 
lows a more quantitative study of these relations and their 
scatter, as well as providing context for our findings on the 
Mass Plane, allowing direct comparison to previous authors. 
All linear ht parameters are given in Table 5. 


4.1 SSP-equivalent parameters versus velocity 
dispersion 

Figure 8 shows the resulting SSP-equivalent parameters for 
our three aperture radii, plotted against the effective velocity 
dispersion, log((Te), given in column 2 of Table 1 in Cappel¬ 
lari et al. (2013a). As shown in the previous section, the pop¬ 
ulation parameters correlate tightly with velocity dispersion, 
conhrming the trends from previous authors. For compari¬ 
son, in Figure 8 we overplot with dashed lines the relations 
from Thomas et al. (2010) based on analysis of 3,360 galax- 
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Figure 8. SSP parameters for the three apertures as a function of log((Te). Magenta points are the galaxies exhibiting low metallicity as 
described in Section 2.3. Dashed lines on the i?e/2 aperture plots show the relations from Thomas et al. (2010), which have a comparable 
aperture size (see text for details). The red line shows the i?e/2 age trend when fitting only galaxies older than 2.5 Gyr, as done by 
Thomas et al. (2010). Solid lines are robust linear fits. Magenta symbols indicate the low-Fe objects highlighted in section 2.5. Solid lines 
are robust linear fits to the black points. Parameters of these fits are given in Table 5. Insert histograms show the residuals about the 
global best fit. 


ies observed with a single fibre by the Sloan Digital Sky 
Survey, selected as part of the ‘MOSES’ project. We make 
this comparison on the Re/"2 aperture, which should be the 
most comparable to the effective aperture of that study. 

The Thomas et al. (2010) lines show general agreement, 
but there are some important differences. In particular, our 
slope with age is significantly steeper, which is trivially ex¬ 
plained by the fact that Thomas et al. (2010) exclude objects 
younger than 2.5 Gyr (at z=0.05-0.06) from their fit. We in¬ 
dicate with a red line our fit applying the same selection, 
and the agreement is excellent. There is also a significant 
difference in slopes for the metallicity relation. Low velocity 
dispersion objects are similar, but the Thomas et al. (2010) 
best-fit linear relation over-estimates the high velocity dis¬ 
persion objects by 0.3 dex - several times larger than the 
total scatter. In Appendix B we show that the difference in 
slope is consistent with the lack of any aperture correction 
applied to the MOSES data to account for the fixed 3” diam¬ 


eter aperture used. This steepens the relation by including 
more of the outer regions in lower-mass (smaller) objects 
relative to the high-mass (larger) objects, which are domi¬ 
nated by the high metallicity inner regions. Even considering 
the offset in slope, the MOSES best-fit linear relation is still 
~ 0.2 dex higher in [Z/H], [a/Fe] also shows a small, less 
significant, offset such that our values are ~ 0.05 dex higher 
than from MOSES. We have verified that these offsets are 
not caused by differences in the models of Schiavon (2007) 
and Thomas et al. (2003, 2004), which, when applied con¬ 
sistently to our line-strengths, show average offsets smaller 
than the la scatter, and which would anyway exacerbate 
the offsets shown in Figure 8. The differences between our 
findings and Thomas et al. (2010) may arise from the differ¬ 
ent combination of indices used by the two studies, but it is 
beyond the scope of this paper to explore this further. 

Contrary to Thomas et al. (2010), we do not find strong 
evidence for bimodality in our distribution of age and veloc- 
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ity dispersion. Our sample size may be too small to constrain 
this, though we note that splitting our sample using their 
age threshold of 2.5 Gyr and our i?e/2 values yields a similar 
fraction of ‘young’ galaxies (23 out of 260, or 8.8%, for our 
sample; 10.15% in Thomas et al. 2010). 

From Figure 8, the effect of different aperture sizes can 
be seen. Most notably, as we ‘zoom out’ from the i?e/8 aper¬ 
ture to Re, the age-CT relation becomes systematically more 
shallow, showing again that young stars, when present, are 
preferentially found in the central regions of the low velocity 
dispersion galaxies. The [Z/HJ-cr relation becomes steeper 
and the average metallicity decreases from around solar for 
Re/S to -0.2 within Re, showing that early-type galaxies are 
only metal rich in their central regions, and are actually 
rather metal-poor objects when considered on the scales of 
one effective radius. Enhancement of alpha elements does not 
show much change considering larger apertures, although we 
remark that the scatter of the [a/FeJ-cr relation increases by 
15% as we go from Re/S to Re, unlike for age and metallicity, 
which become 15% and 10% tighter, respectively. 

4.2 Mass-weighted parameters versus velocity 
dispersion 

Figure 9 presents the mass-weighted age and total metallic¬ 
ity derived from the spectral fitting technique as described 
in Section 2.5. The spectrum fitted is integrated within an 
effective radius, making these values directly comparable to 
the Re aperture of the previous sections. The general trends 
are similar: age and metallicity increase with velocity dis¬ 
persion. However, the correlation with mass-weighted age is 
significantly tighter than with SSP age, showing a factor 2.6 
less scatter (though with a somewhat smaller range in age). 

We also show in Figure 9 the median-binned mass- 
weighted ages to illustrate better the general trend, which 
suggests a possible break in the slope of the relation at 
velocity dispersions higher than ~ 160 kms“^. We caution, 
however, that this is largely a result of limiting the maximum 
age of the stellar population models to 14 Gyr. Including 
the maximum-aged MIUSGAT templates (17.8 Gyr) in the 
star formation history removes any significant deviation 
from a linear relation of mass-weighted age with log(CTe). 
In addition, by fitting a linear relation to points below 
the apparent break, the extrapolated age and scatter at 
the highest velocity dispersion is 18 ± 4 Gyr - only mildly 
inconsistent with the fiducial age of the Universe. This 
suggests that a linear relation is indeed an appropriate 
choice. We return to this issue in Section 4.4. 

The metallicity trend is remarkably similar considering 
the vastly different approaches and models used. These re¬ 
sults give direct empirical confirmation of the fact that early- 
type galaxies are overwhelmingly composed of old stars, but 
in general have a small (but non-negligible) mass-fraction of 
younger stars that contribute disproportionately to the in¬ 
tegrated light of the galaxy - a fact inferred from early SSP 
analyses (Gonzalez 1993; Trager et al. 2000) and shown un¬ 
ambiguously with GALEX UV colours (Kaviraj et al. 2007). 
We consider this more quantitatively in Section 5.1. 

The almost negligible change in the metallicity trend 
also confirms the insensitivity of SSP-equivalent metallicity 
to mixed populations as discussed above. The main differ¬ 
ences are a slight reduction in scatter and a steepening of 
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Figure 9. Mass-weighted age, metallicity and half-mass forma¬ 
tion time of the ATLAS®® sample plotted against velocity disper¬ 
sion. Magenta points are the Fe-poor objects as per Fig. 8. Insert 
histograms show the residuals about the global best fit. Solid lines 
are robust linear fits to the black points. Parameters of these fits 
are given in Table 5. 


the relation, likely due to reduced uncertainties by fitting 
the whole spectrum, and a reduced bias in the contribution 
of young stars responsible for boosting the metallicity in the 
SSP relations. 


4.3 SSP parameters versus dynamical mass 

Figure 10 (left) presents the SSP population trends as mea¬ 
sured against Mjam for the IRe aperture (the smaller aper¬ 
tures follow similar relative differences as those mentioned 
in the previous section). We find that the general trends of 
population parameters are consistent with the picture de¬ 
rived from velocity dispersion alone, such that lower mass 
galaxies tend to show younger ages, lower metallicities, and 
lower abundance ratios. 

The relations with mass are generally more shallow com¬ 
pared to the trends with Ue, especially in age and metallic¬ 
ity, owing to two key differences. Firstly, galaxies with the 
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Figure 10. Left: SSP parameters within one effective radius as a function of dynamical mass, Mjam- Right: Mass-weighted age, metallicity, 
and half-mass formation time, as a function of Mjam- Magenta points indicate the Fe-poor objects as in previous plots. Solid lines are 
robust linear fits to the black points. Parameters of these fits are given in Table 5. Insert histograms show the residuals about the global 
best fit. Green points indicate the median values within bins of 20 neighbouring points, with error bars indicating the ‘error’ on the 
median, computed as ct/a/IV, where cr is the dispersion within the bin, and N = 20. The red line indicates a double power law fit to the 
binned values, as described in Section 4.4. 


lowest masses span almost the full range of ages and metal- 
licities, whereas at the lowest velocity dispersion, galaxies 
are generally the youngest and most metal poor of the sam¬ 
ple. Secondly, the high mass regime where galaxies span a 
narrow range of age and metallicity (Mjam> ac¬ 

counts for about one third of the mass range in the sample. 
For velocity dispersion, this low-scatter region (log(CT) > 2.3) 
is less than a quarter of the sample range. When we move 
to the mass-based relations, the velocity dispersion distribu¬ 
tions are essentially ‘stretched out’ among the older objects, 
reducing the slope of the relations. Both effects are a di¬ 
rect consequence of the distribution of galaxies on the Mass 
Plane. 

Previous studies have found tighter relations of stellar 
population parameters with velocity dispersion than with 
mass (Graves & Faber 2010; Wake et al. 2012), as we 
also qualitatively highlighted on the Mass Plane. Insert his¬ 


tograms in Figures 8 through 10 show the dispersions around 
our linear fits and quantify their standard deviations (see 
Table 5 for the fit parameters and standard deviations). We 
likewise find that the SSP residuals are somewhat tighter 
using CTe than when using mass, although we note that the 
formal difference is small (0.01 dex in age and 0.02 in [Z/H]). 
We note that the improvement of Ue over mass as a linear 
indicator of population parameters is significantly more ap¬ 
parent on the Mass Plane directly than by comparing the 
residuals of it’s edge-on projections, where intrinsic scatter 
makes the interpretation less clear. 

4.4 Mass-weighted parameters versus dynamical 
mass 

Figure 10 (right) presents the mass-weighted age, metallicity, 
and half-mass formation time tso, as a function of Mjam- 


© 0000 RAS, MNRAS 000, 1-33 






















































































15 


We see similar trends as before, but there is now a rather 
clear change in the behaviour of age as a function of mass. 
Rather than being approximated as a linear increase of age 
with mass, albeit with a large scatter at lower masses, the 
trend of mass-weighted age and mass shows a distinct non¬ 
linearity. This is made more apparent by considering the 
median value within mass bins, as indicated by the green 
points in Figure 10. The binned mass-weighted ages show 
a distinct departure from the linear fit, whereas the binned 
SSP values remain reasonably approximated by the best fit 
line. 

To characterise this bimodal trend in more detail, we 
fit a double power law on the median-binned mass-weighted 
values of the form: 


log(Age) = 2^^ log(Agei,) 




(4) 


based on the ‘Nuker’ law for surface brightness profile fitting 
(Lauer et al. 1995). The best fit curve has a low-mass slope 
7 = —0.18, a high-mass slope of /3 = —0.016, and a char¬ 
acteristic mass of Mb = 3.1 x 10^° Mq at a ‘break’ age of 
Ageb = 11.8 Gyr that marks the transition of the two slopes, 
with a ‘sharpness’ parameter for the break in slope given by 
a = 10. This bimodal trend is also evident in the other age 
relations, but to a lesser extent. Using the mass-weighted 
age, the presence of a narrow sequence of old galaxies can 
be seen above the transition mass, with a distinct change 
in the mass-scaling behaviour below this mass. This transi¬ 
tion mass is the same as the inflection point of the ‘zone of 
exclusion’ (ZOE) determined by Cappellari et al. (2013b), 
which defines the low-i?™’^-’ envelope of the Mass Plane as 
shown in Figs. 6 and 7. It also corresponds to the character¬ 
istic mass marking the transition between star-forming and 
passive galaxies found by Kauffmann et al. (2003), as well 
as a possible change in slope of the relation between super- 
massive black hole masses and their host galaxy mass (Scott 
et al. 2013). Considering the population trends on the Mass 
Plane, the correspondence of this ‘knee’ in the age-mass re¬ 
lation to the inflection of the ZOE on the mass-size plane 
becomes obvious: below this mass, the low-i?™**^ envelope 
turns toward larger galaxies, with correspondingly younger 
ages and longer formation timescales. 

We note again that we impose an upper age limit of 
14 Gyr on the model templates used, in accordance with the 
fiducial age of the Universe (13.798 ± 0.037 Gyr; Planck Col¬ 
laboration, Ade et al. 2013). This puts a ‘ceiling’ on the 
mass-weighted age values, as discussed previously in Section 
4.2 for the trend of mass-weighted age with log((7e) (see Fig¬ 
ure 9). However, unlike the trend with log((7e), including the 
full available age range of the MIUSCAT spectral models (up 
to 17.8 Gyr) in the spectral fit does not remove the bimodal¬ 
ity. There is a slightly increased spread in mass-weighted age 
at the oldest values, but the mass at which the ages transi¬ 
tion to a narrow range of old values is unchanged, and the 
clear bimodal shape remains. Moreover, extrapolating the 
linear trend in Figure 10 from masses below the transition 
mass results in strongly unphysical ages (« 30 Gyr) at our 
highest masses, suggesting that the departure from a simple 
linear trend is not artificial. 

Mass-weighted metallicity shows a tight relation with 
Mjam, with only slightly higher scatter than with velocity 
dispersion. This supports the interpretation from the Mass 


Plane that iso-surfaces of constant metallicity run steeply 
with mass, especially for mass-weighted metallicity. 

The half-mass formation time is distributed similarly to 
the mass-weighted age, being strongly related to that prop¬ 
erty. We include it here for completeness. 


4.5 Kinematic class 

One of the key findings of the SAURDN survey was the discov¬ 
ery that ETGs have essentially two categories of global dy¬ 
namics, characterised as ‘fast’ and ‘slow’ rotators (Emsellem 
et al. 2007). Fast rotators have extended regular rotation 
fields well described as inclined disks, though they are clearly 
not ‘thin-disk’ systems, having significant velocity dispersion 
also. Slow rotators do not show regular rotation, and instead 
show counter-rotating structures, decoupled cores, or no ro¬ 
tation at all. In Emsellem et al. (2011), we expand on this 
work with the complete ATLAS^° sample, using the quan¬ 
tity Ar as a way of separating these two classes of ETG. 
In Krajnovic et al. (2011), we have shown that, applying 
the Eourier analysis technique of kinemetry (Krajnovic et 
al. 2006), we can further classify the velocity helds of galax¬ 
ies, giving ‘regular’ rotators (corresponding very closely to 
the fast rotators, which show disk-like rotation), and ‘non¬ 
regular’ rotators (corresponding closely to the slow rotators, 
which show various substructures that cannot be described 
by simple disk-like rotation). 

Figure 11 presents the SSP-Ue relations as before, but 
this time indicating the main kinematic classes described in 
Krajnovic et al. (2011). From this figure it can be seen that 
the various kinematic classes span broad ranges in veloc¬ 
ity dispersion and stellar population properties. Non-regular 
rotators are on average slightly older than regular rotators, 
although we note there are several examples of non-regular 
rotators among the youngest galaxies in our sample. There 
is little discernible difference in abundance ratios, though we 
note that the massive non-rotators are among the most en¬ 
hanced objects in our sample. The most noticeable difference 
between the regular and non-regular rotators is with metal¬ 
licity, where the non-regular rotators appear to fall generally 
on or below the global trend, giving rise to a moderately sig¬ 
nificant difference in the residual distributions around the 
global linear fit as quantified by the K-S test probability, 
Po, indicated in the insert histograms. This implies that 
non-regular rotators have experienced less metallicity enrich¬ 
ment, on average, than their regular-rotating counterparts of 
the same mass. In Young et al. (2011), we also found that 
non-regular rotators typically have lower mass fractions and 
detection rates of molecular gas (though there are notable 
exceptions at low velocity dispersion values, e.g. NGG1222, 
NGG3073). Less molecular gas implies less star formation, 
which is the mechanism for enriching the stellar population 
with metals. 

The offset in metallicity at fixed velocity dispersion can 
be interpreted conversely as an offset in velocity dispersion 
for a fixed metallicity, which could arise as a result of a ma¬ 
jor merger. The metallicity offset is also seen in the relation 
with Mjam. Interpreting the offset as one of increased mass 
and velocity dispersion at fixed metallicity, we find that the 
linear-fit mass- {a-) metallicity relation is offset by 0.3 dex 
in log(MjAM) corresponding to a doubling of the mass, with 
a corresponding shift of 0.1 dex in log((T). These shifts are 
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Figure 11. As per Fig 8, but with symbols indicating kinematic 
classes as defined in Krajnovic et al. (2011). Red circles show lit¬ 
tle or no rotation, green ellipses have complex velocity maps with 
no specific features, green triangles are KDCs, orange lemniscates 
have off-centre dispersion peaks indicative of superposed counter¬ 
rotating disks. These correspond to the respective classes a, b, c, 
and d in Krajnovic et al. (2011), referred to collectively as ‘non¬ 
regular rotators’. Purple spindles are the class e from Krajnovic 
et al. (2011), or ‘regular rotators’. Filled spindles have obvious 
bars. Red and purple lines show fits to only the non-regular and 
regular rotators respectively. The solid black lines show fits to all 
classes together. Parameters of these linear fits are given in Ta¬ 
ble 5. Insert histograms show the residuals of the regular (purple) 
and non-regular (red) rotators about the global best fit to the 
combined population, quoting also Pq: the probability that the 
chance deviation between the two samples should be larger than 
that measured, as derived from a K-S test of the two histograms. 
Small values of Pq imply the distributions are different. The only 
significant difference is in the metallicity relation, where the aver¬ 
age offset of the non-regular rotators is consistent with the mass 
increase expected from a single major (equal mass ratio) merger. 


consistent with the predictions of Hilz et al. (2012) for a sin¬ 
gle 1:1 dry merger, where the velocity dispersion is found to 
increase by a factor 1.2-1.5 following a 1:1 merger. This is 
different to the behaviour expected from the idealised virial 
predictions of Naab et al. (2009), where velocity dispersion 
remains unchanged. The idealised case, however, ignores the 
effects of violent relaxation, which acts to unbind a non- 
negligible fraction of the mass and increase the central den¬ 
sity and velocity dispersion. This interpretation implies that, 
on average, non-regular rotators are consistent with having 
experienced a single dissipationless major merger that could 
have shaped their current kinematic characteristics. This pic¬ 
ture is consistent with Khochfar et al. (2011), where semi- 
analytic models show that slow/non-regular rotators have on 
average experienced more major mergers than fast/regular 
rotators. It is also consistent with Scott et al. (2013), where 
we showed that the radial relation of Mg b line strength and 
escape velocity for our sample limits the number of possible 
major mergers to < 1.5, beyond which objects would become 
inconsistent with the observed trends. 

We add a note of caution to this simple interpretation. 
Firstly, if the offset is in the velocity dispersion, one may 
reasonably expect to see some sign of a corresponding offset 
in the other population parameters, which is not observed. 
This is complicated, however, by the differing response of the 
parameters to mass- and luminosity-weighting. Metallicity 
is essentially a mass-weighted property, suggesting that the 
offset is driven by a process that impacts most of the stellar 
mass in these galaxies. Age, however, is sensitive to small 
amounts of recent star formation, making trends difficult to 
interpret. Figure 4 showed that ages > 8 Gyr are close to the 
mass-weighted values. Considering only galaxies with ages 
older than this, the offset in metallicity remains apparent 
in the K-S test, but the relations with age and abundance 
become sufficiently shallow that any offset in age or velocity 
dispersion is poorly constrained, so we cannot firmly rule out 
a shift in the velocity dispersion axis on this basis. 

Secondly, as we have shown in Naab et al. (2014), the 
net angular momentum of a galaxy is a sensitive tracer of its 
evolution, but with a broad range of possible paths. Some 
of the non-regular rotators are among the very metal-poor 
(Fe-poor) objects, which have high mass fractions of cold 
gas and young ages, consistent with a gas-rich accretion or 
merger event. On an individual galaxy basis, the specific 
assembly history can be varied. However, here we highlight 
the offset in the mass- (ct-) metallicity relation as the only 
significant average difference between the stellar population 
properties of these broad kinematic classes. 


4.6 The role of environment 

The ATLAS^® sample spans a range of environment, from 
sparse fields to the central regions of the Virgo cluster. Dif¬ 
ferences in the galaxy populations of high and low density 
environments are well documented (e.g. Thomas et al. 2005). 
However, reported trends and results are surprisingly weak 
considering the predictions of galaxy formation simulations, 
where dense regions differ in the star-formation properties 
and assembly histories (Dekel & Birnboim 2006; Khochfar 
& Ostriker 2008; Cattaneo et al. 2008). 

Figure 12 shows the SSP-equivalent population parame¬ 
ters and velocity dispersion for our Re aperture, but this time 
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we colour Virgo members red, and non-members blue. In¬ 
sert histograms show the residuals of these two sub-samples 
about the global best fit to the combined sample, quoting 
also Pq: the probability that the chance deviation between 
the two samples should be larger than that measured, as 
derived from a Kolmogorov-Smirnov (K-S) test of the two 
histograms. Small values of Pq imply the distributions are 
different. 

Both Virgo and non-Virgo galaxies show similar trends 
of population properties with velocity dispersion. However 
Virgo objects are older by approximately 2.5 Gyr on aver¬ 
age. The two samples deviate in age most significantly when 
comparing the largest aperture size. This reflects the fact 
that when Virgo galaxies have young stars, they tend to be 
centrally concentrated, whereas outside the cluster environ¬ 
ment, recent star formation can be widespread across the 
galaxy on Re scales. Differences in abundance ratios are also 
statistically most significant within the Re aperture, showing 
a very small probability that the measured deviation of the 
two distributions occurred by chance. Virgo galaxies have 
systematically higher abundance ratios on these large spa¬ 
tial scales than their counterparts in the field at the same 
velocity dispersion. At face value, this indicates a truncated 
star formation history on large scales compared to outside 
the cluster, indicative of a process acting to suppress star for¬ 
mation preferentially on larger scales, such as ram-pressure 
stripping (Gunn & Gott 1972). The suppression of spatially 
extended star formation in Virgo cluster early-type galaxies 
is also supported by the study of current-day star formation 
in early-type galaxies by Shapiro et al. (2010), who found 
that all early-type galaxies with ‘widespread’ star formation 
are found outside of dense environments. Likewise, in pre¬ 
vious papers of the ATLAS®® series we have reported on 
truncated distributions of molecular gas (Davis et al. 2013) 
and a lack of extended HI disks (Serra et al. 2012) within 
Virgo compared to lower-density regions, further supporting 
this picture. 

The metallicity distributions are essentially the same 
between Virgo and non-Virgo objects in all aperture sizes, 
though we note that the ‘low-Fe’ objects are not present in 
Virgo, suggesting that the cluster environment specifically 
prohibits the accretion of low metallicity gas, even though 
the incidence of molecular gas detections within Virgo is sta¬ 
tistically indistinguishable from the field (Young et al. 2011). 
We return to this point when considering the empirical star 
formation histories in Section 5.2 

Our findings here are consistent with our spatial gradi¬ 
ent analysis presented in Scott et al. (2013). Interestingly, al¬ 
though the global properties show differences between Virgo 
and non-Virgo objects, objects in both environments show 
the same radial gradient of Mg 6 with local escape velocity 
(though Virgo is offset to higher Mg 6 values). This implies 
that the mechanism acting to modify the star formation pro¬ 
cess in the cluster galaxies maintains the close relationship 
between the stellar population and the local gravitational 
potential. 

Virgo is the largest bound group of galaxies in our sam¬ 
ple. However, there are regions of our volume that have 
small-scale densities comparable to parts of Virgo. Like¬ 
wise the outer regions of Virgo have densities comparable 
to that in the field. In Cappellari et al. (2011b), we showed 
that using a local estimator of environmental density gives 
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Figure 12. Same as Fig. 8, but separated into Virgo member 
galaxies (red symbols) and non-Virgo (blue). Magenta cross-bars 
show the outlying Fe-poor galaxies. Black, red and blue lines show 
the error-weighted best fit straight line to all galaxies (excluding 
Fe-poor objects), Virgo, and non-Virgo galaxies respectively. In¬ 
sert histograms show the residuals of the best fit line to all galax¬ 
ies, separated into Virgo (red) and non-Virgo (blue) members. 
Pq gives the probability that the two histograms are drawn from 
the same populations, as derived from a K-S test. Metallicities 
are very similar, but age and [«/Fe] show a high probability of 
sampling two different distributions. 


a well-defined relation of kinematic morphologies with envi¬ 
ronment. Here we explore if such similar trends exist for the 
stellar populations. 

Figure 13 shows the LOESS-smoothed distribution of 
SSP-equivalent and mass-weighted population parameters 
(both use Re apertures) as a function of dynamical mass 
and environmental density. Environment is quantified by 
E 3 = Ngea/iT^Rs), which gives the surface density in Mpc“® 
of galaxies inside a cylinder of radius R 3 and height h = 
600 kms“® (i.e. AVhei < 300 kms“®) centered on the galaxy, 
which includes the Vgai = 3 nearest neighbours from the par¬ 
ent sample, as described in Cappellari et al. (2011a). 

It can be seen from Figure 13 that regions of constant 
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Figure 13. Distribution of dynamical mass Mjam versus environment density tracer E 3 , defined in the main text. In each panel, points 
are coloured according to the stellar population parameters indicated in the plot title, after locally averaging the measured values using 
the LOESS technique. The top row provides the SSP-equivalent parameters of age, metallicity and abundance ratio; the bottom row gives 
the mass-weighted age, metallicity and half-mass formation time. The far-right histograms indicate how Virgo cluster members (red) and 
non-members (black) are distributed with E 3 . 


stellar population parameters show varying dependency on 
both mass and environment. The population properties are 
essentially independent of environment over two orders of 
magnitude (within the range —1 ^ log(E 3 ) ^ 1). The 
highest-density regions, dominated by Virgo cluster mem¬ 


bers (see vertical histograms for the distribution of members 
with log(E 3 )), shift markedly to older ages, higher abun¬ 
dance ratios, and earlier half-mass formation timescales for a 
given mass with respect to the lower density regions. Metal¬ 
licity shows only a weak difference between high and low 


© 0000 RAS, MNRAS 000 , 1-33 










































19 


Table 5. Line fit parameters. 





(i) SSP = 

alog((Te)+f) (Figure 8) 




log(Agessp) 
a b 

(7 

a 

[^/-ffjssp 

b 

CT a 

[o/Fe]ssp 

b 

CF 

Notes 


1.14± 0.09 

-1.59± 0.13 

0.22 

0.31± 0.04 

-0.62± 0.06 

0.10 

0.37± 0.03 

-0.61± 0.04 

0.06 

Re/8 

0.93± 0.08 

-1.07± 0.11 

0.20 

0.45± 0.04 

-1.09± 0.06 

0.10 

0.33± 0.03 

-0.49± 0.04 

0.07 

Re 12 

0.83± 0.08 

-0.83± 0.11 

0.19 

0.45± 0.04 

-1.16± 0.06 

0.10 

0.31± 0.03 

-0.44± 0.05 

0.08 

Re 


(ii) SFH = alog((Te)+fe (Figure 9) 


log(AgesFH) 

a b cr 

a b cr 

tso 

a b (j 

Notes 

0.54± 0.04 -0.12± 0.06 0.10 

0.63± 0.04 -1.52± 0.06 0.09 

-10.47± 0.76 25.30± 1.11 1.93 

Re 


(iii) SSP = alog(MjAM)+l> (Figure 10, left) 

log(Agessp) 

a b cr 

[Z/^]SSF [o^/-^e]ssP 

a b cr a b (j 

Notes 

0.25± 0.03 -1.70± 0.10 0.21 

0.12± 0.02 -1.49± 0.05 0.11 0.08± 0.01 -0.66± 0.04 0.08 

Re 


(iv) SFH = alog(ALjAM)+^ (Figure 10, right) 


log(AgesFH) 

a b cr 

a b cr 

a 

tso 

b a 

Notes 

0.15± 0.02 -0.60± 0.05 0.12 

0.21± 0.01 -2.39± 0.05 0.09 

-2.94± 0.32 

34.45± 1.04 2.13 

Re 


(v) SSP = alog((Te)+6 (Figure 11, kinematic classes) 


log(Agessp) 
a b 

cr 

[Z/H]ssp 
a b 

cr 

a 

[a/Fe]ssp 

b 

cr 

Notes 

0.97± 0.09 

-1.12± 0.12 

0.20 

0.53± 0.05 

-1.34± 0.07 

0.10 

0.24± 0.04 

-0.28± 0.06 

0.09 

Regular Rotators 

0.82± 0.14 

-0.81± 0.20 

0.17 

0.74± 0.10 

-1.86± 0.14 

0.12 

0.17± 0.07 

-0.14± 0.10 

0.09 

Non-Regular Rotators 

0.93± 0.07 

-1.06± 0.11 

0.20 

0.57± 0.04 

-1.43± 0.06 

0.12 

0.22± 0.03 

-0.25± 0.05 

0.09 

All (inc. low-Fe) 


(vi) SSP = alog((7e)+fe (Figure 12, Virgo members) 


log(Agessp) 
a b 

cr 

[Z/H]ssp 
a b 

cr 

a 

[a/Fe]ssp 

b 

cr 

Notes 

0.51 ± 0.11 -0.06 ± 0.16 

0.96 ± 0.09 -1.14 ± 0.13 

0.14 

0.18 

0.48 ± 0.07 -1.23 ± 0.10 

0.44 ± 0.05 -1.14 ± 0.07 

0.08 

0.10 

0.25 ± 0.05 
0.36 ± 0.04 

-0.30 ± 0.07 
-0.57 ± 0.06 

0.06 

0.08 

Members 

Non-Members 



Note: The above table gives the parameters for the linear fits shown in the figures indicated in the table. The fits were made using the 
outlier-resistant two-variable linear regression method implemented in the IDL routine ‘robustjinefit.pro’ of the IDL NASA Astronomy 
Library (Landsman 1993). The a column gives the standard deviation of the fit to the ATLAS^^ data. The ‘Notes’ column indicates the 
aperture size and relevant sub-samples used in the fit. If no aperture size is given. Re is assumed. 


densities, in the sense that the highest densities show mildly 
higher metallicities. At the lowest densities (log(E 3 ) < —1), 
there is a weak inflection towards even younger ages, lower 
abundances, and longer half-mass formation timescales. 

Most of our sample (62%) lie within the range — 1 
log(S 3 ) ^ 1, within which galaxy mass is the only significant 


driver of stellar population differences. As we have shown, 
mass-weighted properties are less sensitive to small amounts 
of recent star formation, and therefore trace the build up 
of stars over the entire lifetime of a galaxy. The consistent 
trends between the mass-weighted and SSP-equivalent in¬ 
dicators therefore imply that, on average, the influence of 
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Figure 14. The average star formation history of the corresponding bin in dynamical mass, Mjam- The star formation history is quantified 
here as the relative fraction of stellar mass formed at each epoch, derived from the regularised weighted combination of MIUSCAT SSP 
models fitted to the aperture spectrum. Hatched regions indicate the standard error on the mean at each lookback time, estimated from 
the standard deviation in each mass-bin of age divided by ^Jn, where n is the number of objects in the mass bin as indicated. The upper 


horizontal scale indicates corresponding redshift, assuming f2o = 


environment on the stellar populations of early-type galax¬ 
ies has a long-term influence on the galaxies’ evolution. The 
suppression of only recent star formation in cluster galax¬ 
ies, for example, is not sufficient to explain the trends we 
see, which extend across all masses, and in galaxies where 
significant levels of star formation is long-since over. As we 
show in Section 5.1, there must also have been a more rapid 
build of up massive galaxies in high density environments 
at early times, followed by subsequent preferential suppres¬ 
sion of star formation in those regions. Outside of cluster¬ 
like over-densities, the mass of individual galaxies (or per¬ 
haps more accurately, the parent haloes of today’s galaxies) 
seems to have been the dominant influence governing the 
star-formation history over the effects of the local galaxy 
density. This supports a picture of separable effects of mass 
and environment on the stellar populations of ETGs, con¬ 
sistent with Peng et al. (2010). The weak influence of envi¬ 
ronment (as compared with mass) in low- and intermediate- 
density regimes is also similar to the hndings of Thomas et 
al. (2010). 


5 STAR FORMATION HISTORIES 

In this section, we use our spectral fits to the integrated Re 
aperture spectra to investigate the history of star formation 


-|- Oa = 1, where 12m = 0.27 and Ha = 0.73. 


within our sample galaxies. The results which follow use the 
same spectral fits used to determine the mass-weighted pop¬ 
ulation parameters presented above. As described in Section 
2.5, we make a regularized linear ht to our aperture spectra 
of SSP model spectra sampling a grid of age and metallic- 
ity. The resulting distribution of weights provides an esti¬ 
mate of the star formation history (SFH) of the galaxy in 
question, by describing the mass contributions of different 
populations with given formation epochs. The spectral fit 
is made using model spectra logarithmically-spaced in age. 
Here we rebin the resulting logarithmically-spaced weights 
onto a linear time axis, rigorously preserving the total mass 
fractions within each linear time step (set as 1 Gyr). 

5.1 Star formation histories as a function of mass 

Figure 14 shows the measured SFHs of our sample aver¬ 
aged within five bins of dynamical mass (as measured at 
z = 0), spanning the range of our sample. The SFH is pre¬ 
sented as the mean fraction of stellar mass formed in each 
1 Gyr step, down to mass fractions of 1%, for galaxies within 
the (present-day) mass bin. The variation of mass fraction 
formed within each 1 Gyr bin is indicated by the hatched 
regions, quantified as the mean standard deviation within 
each 1 Gyr bin. 

On average, objects with Mjam> IO^^Mq (> 90% of 
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Figure 15. Cumulative functions of the star formation histories shown in Figure 14, showing the integrated mass fraction existing at 
given time. Colours correspond to mass bins as in Figure 14. Hatched regions indicate the dispersion of values at each time step. 


our sample) had their peak in star formation at the earli¬ 
est epochs, followed by a decline that is systematically more 
rapid at higher (present-day) masses. In our lowest mass 
bin, the behavionr is different: on average, these objects 
build a gradually increasing fraction of their stellar mass 
as the universe evolves, until a more rapid decline beginning 
around 4 Gyr ago. This decline is most easily interpreted 
as a result of our sample selection, which targeted galaxies 
consistent with an early-type morphology, and thus prefer¬ 
entially selects objects with low present-day star formation 
rates (though this is not a selection requirement, and many 
objects show evidence for ongoing star formation). We note 
also that this lowest mass bin likely suffers from incomplete¬ 
ness, as the sample was selected based on K-band luminosity, 
which has some variation in mass-to-light ratio with stellar 
population properties. This may bias against old populations 
at this low mass extreme, as they will be slightly fainter than 
their younger counterparts. Nonetheless, the lowest mass bin 
is consistent with the overall trend of more extended star for¬ 
mation in lower mass early-type galaxies, with a truncation 
at more recent times than the high mass objects. 

The star-formation histories shown in Fig. 14 provide 
an alternative picture to that of the schematic view inferred 
from assuming linear trends of SSP-equivalent parameters, 
as in Thomas et al. (2005, 2010). In that case, galaxies of a 
given mass have star formation rates that rise and fall sym¬ 
metrically in time, being described by a Gaussian centred 
on the mean SSP-equivalent age, with a width given by a 
conversion of abundance ratio to star-formation timescale. 


With that approach, the stellar population scaling relations 
arise from a systematic delay in the peak of star formation 
in lower mass galaxies towards later times, combined with a 
corresponding increase in the star-formation timescale. 

By comparison, our non-parametric approach shows 
that objects of all masses experience a significant epoch 
of star-formation at the earliest times, and that the subse¬ 
quent variation in present-day stellar population properties 
is mainly driven by the differing rate at which star formation 
declines for galaxies of different masses. The overall picture is 
qualitatively similar between both approaches. However, by 
reproducing the observed spectra with a distribution of ages 
and metallicities directly, our non-parametric star-formation 
histories emphasise the significance of early star formation 
in all galaxies, and the subsequent decline in star formation 
activity that varies systematically with mass, independent of 
the timescale of the individual star formation events. This re¬ 
flects the conceptual differences of fitting the full spectrum 
with a distribution of models to infer the star formation 
history directly, compared to using line-strength indices to 
estimate mean properties. Our findings are in good agree¬ 
ment, for example, with the previous studies of Heavens et 
al. (2004) and Banter et al. (2007), who used comparable 
spectral-fitting methods on thousands of aperture spectra 
from SDSS (though without any morphological selection ap¬ 
plied). 

Figure 15 presents the associated mean cumulative mass 
fraction as a function of look-back time for the same mass 
bins as in Figure 14, as well as the total cumulative formed 
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mass for the entire sample. Assuming that our sample rep¬ 
resents most of the stellar content of today’s local Universe, 
approximately 50% of all stars formed within the first 2 Gyr 
following the big bang. The vast majority of these stars re¬ 
side today in the most massive galaxies, which themselves 
formed 90% of their stars within the first 3 Gyr (i.e. by « ~2). 
The low mass objects, in contrast, have formed barely half 
their stars in this time interval, forming the remaining 50% 
of their stellar mass in the last 6-10 Gyr. We note that, from 
this analysis, we can only tell when the stars formed, and 
not whether they formed within the gravitational potential 
of the current-day galaxy, or instead were formed in sepa¬ 
rate, smaller systems that later combined to form the galaxy 
we see today. 

Figure 16 shows the fractional formation masses in ab¬ 
solute terms (top panel), which allows comparison to other 
studies that measure star formation rates. We do this by 
assuming that our star formation histories apply to the 
galaxy’s total mass as quantified by Mjam, scaling this 
present-day mass by the weights returned by pPXF. We also 
derive the history of specific star formation rates (middle 
panel) by normalising the star formation rate by the inte¬ 
grated mass formed by the epoch being considered (though 
the galaxy may not have assembled yet). In addition, since 
our sample is volume-limited, we also present the history of 
star formation rate density for our local volume implied by 
our analysis, both as a function of mass, and the combined 
total for our full sample (bottom panel). 

The top panel of Figure 16 clearly shows that the most 
massive current-day galaxies had the highest absolute star 
formation rates of all such objects for the first 4 Gyr follow¬ 
ing the big bang. The bottom panel of Figure 16, however, 
shows that these objects do not dominate the history of in¬ 
tegrated star formation rate density, due to their smaller 
number density. 

The middle panel of Figure 16 presents the time- 
averaged specific star formation rate (sSFR), defined as the 
mean star formation rate during the 1 Gyr time bin divided 
by the mass formed up to that epoch (taken as the mass 
formed up to the centre of the time bin). For comparison, we 
show direct measurements of the instantaneous sSFR taken 
from the literature. Our derived sSFR values lie generally 
well below the direct measurements, which is most easily 
explained by the inconsistent timescales probed by the two 
measures of star formation rate. Our ‘archeological’ values 
probe the mean time-averaged star formation during each 
10® year time step, whereas the direct measurements probe 
the star formation occurring on timescales of 10® —10® years. 
The super-solar abundance ratios of our most massive (and 
oldest) galaxies require star formation to have occurred on 
characteristic timescales of order 10® year (Thomas et al. 
2005), which is below the time resolution of our spectral fit¬ 
ting at these early epochs. Such shorter ‘bursts’ of star for¬ 
mation would boost our effective star SFR and sSFR values 
by an order of magnitude, potentially making them more 
consistent with direct instantaneous measurements. How¬ 
ever, we cannot directly constrain this timescale during the 
oldest age bins using our current approach. 

In addition to a redshift dependence, the sSFR histories 
show a strong mass-dependence of the sSFR at a given red- 
shift, irrespective of the assumed star-formation timescale. 
This indicates the influence of a mass-dependent feedback 


mechanism that is relatively more eflicient in higher-mass 
galaxies, reducing the sSFR more rapidly than at lower 
masses, and bringing higher-mass galaxies off of the star- 
formation main sequence at earlier times. Through our 
‘archeological’ approach, we can directly trace when galax¬ 
ies turn off of the star-formation main sequence and be¬ 
come quiescent. Due to the problems of comparing instanta¬ 
neous and archeological sSFR estimates, however, we cannot 
easily compare with the measured sSFR main sequence at 
high redshift. Instead, we adopt the sSFR of 0.1 Gyr“^ mea¬ 
sured in local galaxies from Peng et al. (2010) as a robust 
threshold between ‘star forming’ and ‘quiescent’ states, and 
measure the epoch at which our galaxies cross this thresh¬ 
old. Figure 17 presents this ‘quenching redshift’ for the dif¬ 
ferent mass bins. We find a near-linear dependence with 
log(MjAM), with galaxies of mass above 10^^ Mq quenching 
their star-formation at redshift 2 = 2 — 4, whereas galax¬ 
ies below 10^® Mq quench below redshift 1, with our lowest 
mass galaxies quenching below 2 ~ 0.2 (in the last IGyr), 
on average. 

The cosmic SFR density for ATLAS®® given in the bot¬ 
tom panel of Figure 16 also shows measurements from Panter 
et al. (2007) derived using a similar technique to ours using 
SDSS spectra of ~ 3 x 10® galaxies, and from the compila¬ 
tion of Gunawardhana et al. (2013) using various measures 
of the (instantaneous) star formation rate. Our estimate of 
the total star formation rate density is significantly lower 
than these other studies, differing by an order of magnitude 
at most recent times, but with better agreement at the earli¬ 
est epochs. We interpret this as a result of our morphological 
selection, which effectively excludes galaxies with strong to¬ 
tal star formation rates at the current time, and which have 
likely contributed significantly to the integrated star forma¬ 
tion budget of the universe over its lifetime. In contrast, the 
literature studies either impose no particular mass or mor¬ 
phological selection (as in the case of Panter et al. 2007), or 
preferentially select galaxies with detectable star formation 
(as in the data compiled by Gunawardhana et al. 2013), both 
of which will act to increase the inferred integrated star for¬ 
mation rate. Gosmic variance may also play a role given our 
modest survey volume, which for example lacks extremely 
massive clusters like Goma. In general, however, we see that 
early-type galaxies have experienced a similarly-paced de¬ 
cline in star formation rate as the total population, though 
with a more rapid drop off and reduction in contributions 
from galaxies of higher masses. 

5.2 Star formation histories as a function of 
environment 

In Section 4.6 we presented evidence for differing stellar pop¬ 
ulation properties in the different environments present in 
our sample. Here we investigate whether this is also recov¬ 
ered directly in the SFHs determined from spectral fitting. 
Figure 18 shows the mean cumulative fraction of stellar mass 
formed as a function of time for members of the Virgo clus¬ 
ter and the rest of the sample, binned by present-day mass. 
The population outside of Virgo (solid lines) show consis¬ 
tently more extended star formation than the objects that 
are Virgo members (hatched regions), with the exception 
of the lowest mass bin, where the extended star formation is 
similar in both populations. This low-mass bin has just three 
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Figure 16. Top: Star formation rates averaged within the indicated (present-day) mass bins. These are calculated by multiplying the mass 
fractions in Figure 14 by the dynamical mass Mjam- Middle: Mean (time-averaged) specific star formation rate in each AtsFH = 1 Gyr 
time step (thick solid lines), obtained by normalising the star formation rate by the mass formed up to the epoch considered. Direct 
measurements of the instantaneous sSFR from the literature are also shown, either as discrete values (black data points) or empirically- 
constrained continuous relations (broken lines). Bottom: Integrated star formation rate density for our local volume, as a function of 
present day galaxy mass, and the integrated total for our sample (black line). Grey points show instantaneous star formation rate density 
measurements from the compilation of Gunawardhana et al. (2013); red points show the measurements of Panter et al. (2007) using a 
similar approach to our analysis. In all panels, hatched regions show the standard deviation of properties within each age and mass bin. 
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Epoch where sSFR<0.1 Gyr 



Figure 17. Redshift at which the sSFR in different bins of galaxy 
mass, MjaMi crosses the value of 0.1 Gyr“^, corresponding to the 
local sSFR of main-sequence star-forming galaxies from Peng et 
al. (2010). The quenching redshift follows an approximately linear 
relation with mass, as indicated by fit (the function is also given). 
Errors are computed using the spread within each mass bin, as 
indicated with hatched regions in Figure 16. 


Our star formation histories suggest that this is not re¬ 
stricted to a recent suppression of star formation. On av¬ 
erage, cluster objects have always exhibited shorter forma¬ 
tion timescales than field galaxies of comparable present-day 
mass. This finding is in good agreement with the predictions 
of Khochfar & Ostriker (2008), where a combination of ram 
pressure stripping and gravitational heating act to reduce 
the duration of star formation in dense environments. The 
corollary to this is that, at early times when cluster galaxies 
were rapidly forming the bulk of their stars, their effective 
star formation rate would have been higher than that of 
galaxies in the field that would eventually match their mass 
in the present day. Such a picture is consistent with direct 
observations at high redshift, where an apparent ‘reversal’ 
of the low redshift behaviour (where star formation is lower 
in high density regions) occurs, and dense environments ex¬ 
hibit enhanced star formation compared to the field (Elbaz 
et al. 2007; Tran et al. 2010), potentially driven by cold flows 
(Khochfar & Silk 2009). This phase would be largely com¬ 
pleted by redshift z ~ 2, by which time the most massive 
galaxies have formed their stars, and galaxies of lower mass 
and in less dense environments become the primary sources 
of star formation for the remaining 10 Gyr. 


Virgo objects, all of which are located in the outer regions 
of the cluster, and are presumably entering Virgo for the 
first time. Overall, these trends directly show the impact of 
environment on the mass growth and star-formation histo¬ 
ries of early-type galaxies, in addition to the mass-dependent 
effects illustrated in the previous section. 

The apparent suppression of extended star formation 
within Virgo compared to objects of the same mass outside 
of the cluster environment is qualitatively consistent with an 
extension of what is found in Virgo late-type galaxies, where 
the intra-cluster medium appears to truncate extended neu¬ 
tral gas disks with subsequent implications on star-formation 
activity (Koopmann & Kenney 2004; Growl & Kenney 2008; 
Ghung et al. 2009). Indeed, even though the detection rate of 
GO in our sample is not statistically different between Virgo 
members and non-members (Young et al. 2011), we have 
shown in Davis et al. (2013) that this present-day reservoir 
of molecular gas responsible for fuelling current star forma¬ 
tion is spatially less extended in the Virgo galaxies, and often 
shows signs of recent interaction and/or truncation. 

Further evidence comes from the kinematic alignment of 
gas in our sample. In Davis et al. (2011) we showed that gas 
in various phases (ionised, molecular or neutral) is almost 
always aligned with the stars in Virgo galaxies, in contrast 
to galaxies in low-density environments where more than 
half show misalignments. This strongly suggests that the gas 
supply in the cluster is limited to internal recycling. Without 
any additional external accretion of cold gas to maintain or 
periodically boost levels of star formation, galaxies in the 
cluster have less extended star formation histories compared 
to those of similar mass outside the cluster environment. 
In less dense environments, accretion of cold gas can still 
occur, as indicated by the high incidence of gas components 
misaligned to the stellar motions. Recall also that none of the 
extremely low metallicity, gas rich objects shown in Figure 5 
are found in Virgo. 


6 DISCUSSION AND CONCLUSIONS 

Our integral-field spectroscopy provides a unique opportu¬ 
nity to study the mass-scaling properties of stellar popula¬ 
tions in early-type galaxies with high accuracy, by virtue of 
minimising the effects of aperture size and slit losses. For 
a large range of galaxy masses and apparent sizes, we im¬ 
plicitly include the effects of line-strength gradients in our 
measurements of population properties by matching scaled 
apertures with effectively no slit loss effects. Our study also 
uses CTe: the velocity dispersion integrated within , which is 
empirically close to the true second velocity moment (wfos)oo 
that appears in the scalar virial theorem (Binney & Tremaine 
2008), and is directly proportional to mass. Moreover, the 
IFU data allows the dynamical mass to be tightly con¬ 
strained by fitting general dynamical models that account for 
non-homology and dark matter. Our sample selection is also 
based on a classical prescription of determining galaxy type, 
based on the presence of spiral arms and/or a significant dust 
disk. In this way, we do not exclude galaxies that have on¬ 
going star-formation on anything other than morphological 
grounds. Finally, we have a large collection of high quality 
measurements on individual galaxies, allowing the intrinsic 
scatter of our relations to be accurately assessed. 

From consideration of the classic Fundamental Plane, 
in Cappellari et al. (2013a) we showed that the intrinsic 
scatter and tilt of this relation is driven by stellar popu¬ 
lation variations, including the stellar initial mass function, 
and that accounting for these effects in the stellar mass- 
to-light ratio brings galaxy masses into very close agree¬ 
ment with the scalar virial relation Mjamoc Here 

we explored the distribution of stellar population param¬ 
eters on the thin plane of {Mjam, - or the ‘Mass 

Plane’ - and show that there are indeed smooth variations of 
age, metallicity and abundance ratio within this plane. By 
consideration of empirical star formation histories derived 
from spectral fitting, we confirm that this reflects trends in 
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Figure 18. Cumulative star formation history for Virgo members (thick hatched regions) compared to outside of Virgo (solid lines) for 
galaxies averaged in bins of dynamical mass. Line colours show the range of dynamical masses considered, as indicated. 


the mass-weighted properties and star formation timescales 
across the Mass Plane. At fixed mass, galaxies with smaller 
sizes are on average older, more metal rich, and more alpha- 
enhanced than their larger counterparts, showing that they 
formed earlier and in a shorter time. These trends tie di¬ 
rectly to a higher mass normalisation of the IMF in these 
more compact objects (Cappellari et al. 2012; McDermid et 
al. 2014), as well as a larger bulge fraction and lower molec¬ 
ular gas fractions (Cappellari et al. 2013b). 

Why are today’s massive galaxies older, more metal rich, 
and more enhanced in alpha elements compared to their 
lower-mass counterparts? There seem to be two processes 
at work. All of today’s early-type galaxies share an early 
period of intense star formation, making the systems com¬ 
pact, metal rich and alpha enhanced. This is consistent with 
the picture of Khochfar & Silk (2006a), who argue that high 
levels of dissipation at high redshift leads to more compact 
galaxies, with consequently older ages and enhanced abnn- 
dance ratios. Lower-mass systems can continue to accrete 
gas from the inter-galactic medium, thereby increasing their 
apparent size and forming new stars, resulting in younger 
average ages when viewed today, with lower average metal- 
licities. Dissipation of accreted gas results in a extended and 
regular stellar rotation field in these objects, which remains 
in tact in the absence of any major merger events, produc¬ 
ing a ‘regular’ or ‘fast’ rotator (Khochfar et al. 2011; Naab 
et al. 2014). In the case where a major merger does occur, 
depending on the orbital parameters of the interaction, the 


resulting system can remain as a ‘fast’ rotator, or in the 
case of misaligned spin vectors of the progenitors, give rise 
to a ‘two-sigma’ galaxy with a significant fraction of stars on 
counter-rotating orbits (Bois et al. 2011). 

Massive systems, in contrast, quench early, and there¬ 
after predominantly accrete stars (Khochfar & Silk 2006; 
Oser et al. 2010). Most of these stars have also formed early, 
but in smaller systems with lower metallicity. Via this stellar 
accretion, these galaxies will also increase in size, but remain 
old on average. The accreted low mass systems are dispersed 
at a radius proportional to where the local potential ap¬ 
proaches the binding energy of the accreted system (Villum- 
sen 1983), resulting in a gradient of metallicity, decreasing 
to large radii (Hirschmann et al. 2013; Navarro-Gonzalez et 
al. 2013). Moreover, the collisionless accretion also reduces 
the average angular momentum of the system, resulting in 
massive ‘non-rotating’ galaxies (Naab et al. 2014). Where 
neither accretion of gas nor stars is efficient, galaxies remain 
compact, with characteristic stellar populations that are rel¬ 
atively old, metal-rich and alpha-enhanced for their mass. 
We note that in both cases (gas accretion at lower mass, 
and stellar accretion at higher mass), the stellar metallicity 
must remain closely related to the local potential, as traced 
by the spatially-resolved relationship between stellar popu¬ 
lation parameters and the local escape velocity (Scott et al. 
2013). 

Cappellari et al. (2013b) discuss evolution in the Mass 
Plane at length, describing the processes by which objects 
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may evolve within that plane. In summary, the evolutionary 
picture described there is one of increasing bulge fraction (or 
decreasing disk fraction) as mass increases, bringing galaxies 
to higher and more compact sizes. This process may be 
driven by dissipative mergers, which encourage central star 
formation and increase the effective concentration of stel¬ 
lar mass. Above a characteristic mass of log(MjAM) = 11.2, 
the range of early-type galaxy sizes is dramatically reduced, 
together with the scatter in stellar population properties. 
The process driving mass growth in this regime is effective 
at increasing galaxy size without signihcantly changing the 
stellar population properties or velocity dispersion. Simula¬ 
tions of collisionless (or ‘dry’) mergers show behaviour that 
is largely consistent with this picture (Khochfar & Burkert 
2003; Naab et al. 2009; Hilz et al. 2013), and given the low 
gas fractions associated with the oldest and/or most massive 
early-type galaxies (Young et al. 2011), are a suitable mech¬ 
anism for driving the mass growth along the upper mass 
envelope of the mass-size plane. Moreover, we present evi¬ 
dence that the mean scaling relations for the metallicity in 
non-regularly rotating galaxies is offset to higher mass and 
velocity dispersion for a given metallicity, consistent with a 
major merger origin for these objects. 

We also derived empirical, non-parameteric star forma¬ 
tion histories (SFHs) for our sample, in order to derive mass- 
weighted mean ages and metallicities. As expected, the mass- 
weighted ages are significantly older than the SSP-equivalent 
ages, which are strongly biased by the presence of bright 
young stars in galaxies with recent star formation. 

In addition to providing mass-weighted (instead of SSP- 
equivalent) stellar population properties, our derived star 
formation histories themselves provide unique insights into 
the formation of our sample galaxies. We hnd a very regular 
dependence of star formation history on present-day mass, 
consistent with similar results using the SDSS (Heavens et 
al. 2004; Panter et al. 2007) and at higher redshift (Juneau 
et al. 2005). Within bins of galaxy mass, the peak of star 
formation is, on average, within the first billion years fol¬ 
lowing the big bang for all galaxies with Mjam> 

This provides an alternative view to the picture of staggered 
Gaussian star formation that is implied assuming linear re¬ 
lations of SSP-equivalent parameters and abundance-ratios 
to derive formation timescales, as in Thomas et al. (2005, 
2010), though the qualitative picture is similar. 

Through consideration of the mean absolute star for¬ 
mation rate, we demonstrate that the classic ‘downsizing’ 
scenario of Cowie et al. (1996) applies to early-type galax¬ 
ies, showing that less massive galaxies have more extended 
star formation histories, and as a result, the characteristic 
mass of the objects dominating the star formation rate de¬ 
creases as the Universe evolves. This trend is driven by the 
increasingly rapid fall off in star formation for galaxies with 
higher masses. 

Our most massive galaxies (11.5 < log(MjAM)< 12) 
have formed 50% of their stars no more than 2 Gyr after the 
big bang (a ~ 3). By comparison, the lowest-mass galaxies 
(9.5 < log(MjAM)< 10) reach the same fraction after around 
8 Gyr {z ~ 0.6). This is in reasonable agreement with recent 
simulation results (e.g. Oser et al. 2010), where models are 
indeed characterized by early intense star formation fed by 
cold flows and subsequent assembly of the formed stars, with 
ongoing ‘in situ’ star formation more prevalent in lower-mass 


galaxies. Although we cannot constrain where the stars have 
formed from our integrated measurements, it is certainly the 
case that the recent (and even ongoing: Shapiro et al. 2010; 
Davis et al. 2014) star formation is more present in the cen¬ 
tral regions of our lower-mass galaxies, supporting the view 
that this is ‘in situ’. 

Finally, we explore the environmental dependence of our 
star formation histories, and find that, on average, galaxies 
in Virgo have an additional suppression of extended-duration 
star formation with respect to galaxies outside the clus¬ 
ter environment. As a result, comparing galaxies of similar 
mass, Virgo members build their stellar mass on a shorter 
timescale than non-members - an effect that is most pro¬ 
nounced at intermediate masses (10^° <Mjam< IO^^Mq). 
The highest mass (Mjam> IO^^'^Mq) galaxies are very alike 
in their star formation histories, forming stars on compa¬ 
rably short timescales (90% of stellar mass formed within 
3Gyr). This suggests that the star formation properties of 
these galaxies are relatively unaffected by their current-day 
environment. Interestingly, galaxies with the lowest masses 
(Mjam< IO^^Mq) also share very similar star formation his¬ 
tories, with protracted formation of stars until recent times. 
These lowest-mass Virgo objects are at large cluster-centric 
distances, suggesting they have only recently entered the 
dense cluster environment for the first time, and the effec¬ 
tive suppression of star formation by the environment has 
not yet made an impact. 
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APPENDIX A: LINE-STRENGTH 
CALIBRATION AND ERRORS FROM 
STANDARD STARS 

For observations made before the installation of the VPH 
grating in SAURQN in 2004, we refer the reader to Kuntschner 
et al. (2006) for details. An extensive set of line-strength 
standard star observations were also obtained post-VPH in¬ 
stallation to provide a calibration to the Lick system of in¬ 
dex measurements, and to monitor run-to-run and inter-run 
variations. The stars were chosen to overlap with several 
well-known stellar libraries, including Worthey et al. (1994), 
Jones & Worthey (1995) and MILES (Sanchez-Blazquez et 
al. 2006). In total, 67 stars were observed, with several re¬ 
peat observations, giving a total number of 160 stellar spec¬ 
tra. Wherever possible, velocity measurements for the stars 
were derived from the matching MILES library spectrum, 
otherwise the closest spectral type from the MILES library 
was used. This minimized template mismatch effects from 
introducing large velocity errors, and subsequent errors in 
the index measurements. The stars do not suffer from the 
field coverage issue mentioned previously, so we include it 
here for reference. 

Figure Al presents the difference between our measure¬ 
ments and the standard reference values of the Lick/IDS 
system (Worthey et al. 1994). Since some stars have several 
repeat measurements, and others do not, we use an error- 
weighted mean to derive the average offsets. Six stars were 
observed four times or more, upon which we base estimates 
of our intrinsic measurement errors in Figure Al and given in 
Table A2. This characteristic error was adopted for each indi¬ 
vidual measurement, and repeat measurements are averaged. 
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Figure Al. Comparison of Lick indices measured from the 67 
standard stars overlapping between the ATLAS®^ survey and 
those from the Lick/IDS system (Worthey et al. 1994). The de¬ 
rived biweight-mean offsets (shown as a dashed line) and errors on 
the mean are given in Table Al. Dots show individual measure¬ 
ments; diamonds show individual and (where available) combined 
measurements. A representative error bar is plotted in the upper 
left of each panel showing the empirical scatter in repeat mea¬ 
surements from six stars having four or more observations over 
the four observing runs. This value is also tabulated in Table Al, 
and shows that the scatter in this comparison is largely driven by 
scatter in the original Lick measurements, which have an average 
uncertainty around twice as large. 


with the characteristic error scaled by (n — where n 

is the number of repeat measurements. We effectively ignore 
the intrinsic errors of the reference values themselves, which 
in fact are the dominant source of scatter (compare the er¬ 
rors in Table A2 with the scatter in Figure Al, quantihed by 
CT in the panel titles). As with all figures in this section, we 
employ outlier-resistant ‘robust’ estimators of the standard 
deviation. 

The average offsets from our post-VPH standard stars 
are consistent within the la uncertainties of the pre-VPH 
‘stars only’ values of Kuntschner et al. (2006, column 2 of 
their table 3), with the exception of Mg b, for which we derive 
a larger offset, indicating slightly weaker index values with 
the VPH grating. To further compare the pre- and post-VPH 
consistency of our SAURON data, we compare measurements 
of the same stars measured with the two different gratings. 
Figure A2 shows this comparison, where the values have 
been averaged wherever possible. The mean offsets in this 
comparison are fully consistent with the differences found 
in the comparison with the Lick/IDS values, indicating that 
the small differences between the two gratings are real. We 
therefore apply different offsets for the pre- and post-VPH 
data. The final applied offsets are given in Table Al. We also 
derive an estimate of our uncertainties from this differential 
comparison, assuming that the scatter comes from equally- 
sized errors from both observations. The inferred values are 
also given in Table A2, and are in good agreement with the 
estimates from repeated stars. 

Whilst these error estimates are inferred from observa- 
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Figure A2. Comparison of Lick indices measured from the 43 
stars observed with SAURON before (from Kuntschner et al. 2006) 
and after the installation of the VPH grating. Repeat measure¬ 
ments were averaged in both data sets, and corresponding errors 
combined. The plotted errors assume the errors come equally from 
both data sets, giving croverlap in Table A2. Dashed line indicates 
the biweight mean. 


Table Al. Estimate of Lick offsets derived from all stars. 


Index 

Pre-VPH Offset 

Post-VPH Offset 

H/3 

-0.13 ± 0.05 A 

-0.09 ± 0.03 A 

FeSOlS 

-1-0.28 ± 0.05 A 

-10.37 ± 0.06 A 

Mgfe 

-1-0.00 ± 0.05 A 

4-0.07 ± 0.03 A 

Fe5270 

-1-0.00 ± 0.05 A 

4-0.07 ± 0.02 A 


Note: For each index, the biweight mean and dispersion estimate 
was used to reduce the influence of outlying values. The error 
on the mean offset derived from the N observations is given as 
the dispersion scaled by l/v/V. Pre-VPH offsets are taken from 
(Kuntschner et al. 2006). 

tions spanning several years, they are still ’internal’ error es¬ 
timates, coming from the same instrument and similar data 
reduction etc. As a way to check our ’external’ errors, i.e. 
how we may compare with others, we present in Figure A3 
a comparison of the Lick indices derived from our measure¬ 
ments (post-VPH) and from the MILES spectra for the same 
stars. The scatter is consistent with our internal comparison, 
suggesting that our error estimates are reasonable. Moreover, 
the mean offset and scatter are comparable to those found 
between targeted stellar libraries (e.g. see Sanchez-Blazquez 
et al. 2006). 


APPENDIX B: COMPARISON WITH OTHER 
AUTHORS 

B1 Comparison of line-strengths 

The dominant uncertainty in line strength measurements 
from high signal-to-noise ratio data such as is available for 
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Figure A3. Comparison of Lick indices measured from the 57 
standard stars in common with the MILES library of Sanchez- 
Blazquez et al. (2006). Small dots show all 147 individual mea¬ 
surements, and diamonds show individual and averaged values 
where repeat measurements are available. The mean, error on the 
mean and dispersion are given in the plot titles, which are com¬ 
puted from the individual values to provide a measure of our un¬ 
certainties. The representative error bar (filled symbol) assumes 
the scatter comes equally from both data sets. 


Table A2. Error estimates derived from the standard star mea¬ 
surements. cTp^epeat is derived from the mean standard deviation 
of six stars with four or more measurements, croveriap is derived 
by assuming constant and equal errors between pre- and post- 
VPH observations of stars in common to both set-ups. ctj^ules 
is derived by comparing our Lick indices with those from MILES 
spectra for the same stars, and assuming the scatter comes from 
both data sets equally. The last column gives the average devia¬ 
tion, which we take as our minimum uncertainty in each index. 


Index 

^Repeat 

'^Overlap 

e'MILES 

(7 

H/3 

0.10 A 

0.08 A 

0.11 A 

0.10 A 

Fe5015 

0.14 A 

0.21 A 

0.19 A 

0.18 A 

Mgb 

0.12 A 

0.08 A 

0.08 A 

0.09 A 

Fe5270 

0.11 A 

0.11 A 

0.11 A 

0.11 A 


nearby galaxies are the systematic errors from data calibra¬ 
tion, rather than statistical errors from noise. We therefore 
choose from the literature two studies that have reasonable 
overlap in sample galaxies, namely the extensive Lick/IDS 
sample of Trager et al. (1998) (90 galaxies in common), and 
the higher-quality data from Gonzalez (1993) (20 galaxies in 
common). 

For the Trager et al. (1998) data, as per Kuntschner 
et al. (2006) we extract a circular aperture with the same 
area as their 1.4''x4"(i.e. with a radius of 1.335"), since 
we do not know the position angle of the slits used. The 
Trager et al. (1998) data is also not corrected for emission, 
which is a major difference to our data. We therefore mea¬ 
sure the comparison values with no correction for emission. 
Figure B1 presents the comparison, with robust mean, er¬ 


0.08+/- 0.03 (a= 0.25) 



Hbeto Troger 



FeSOI 5 Troger 

0.09+/- 0.03 {a= 0.29) 



Figure Bl. Comparison with the Lick indices of Trager et al. 
(1998, T98). No correction for emission is made in either data set. 
The T98 data were measured in 1.4" x 4" apertures of unknown 
orientation. The ATLAS^^ data (A3D) are measured in circles of 
equivalent area (i.e. of radius 1.335”). Errors show the quadrature 
sum of both data set errors. Plot titles show the robust mean (also 
shown as a dashed line), error on the mean, and robust dispersion 
(in parentheses). 


ror on the mean, and dispersion (in parenthesis) shown in 
the title. Observational errors for both data sets have been 
added in quadrature. The comparison is reasonable, though 
systematic differences can be seen, in particular an apparent 
trend with increasing Fe5015. Corrections for velocity disper¬ 
sion are comparable between the two studies to ~ 1%. The 
outlier in the H/3 comparison, NGC3032, demonstrates the 
sensitivity to such comparisons to PSF and aperture differ¬ 
ences. This particular galaxy has an unresolved blue nucleus 
where H/? strengthens considerably (McDermid et al. 2006), 
making the PSF a dominant uncertainty. 

Figure B2 presents the comparison with the i?e/8 aper¬ 
ture values from Gonzalez (1993), with titles and errors as 
per Figure Bl. In this work, careful consideration of the aper¬ 
ture definitions was made, and the translation between the 
slit radial length and the corresponding circular aperture 
size are conveniently tabulated there (his table 4.7). Using 
this radius, we extract spectra within circles, and compare 
the measured line-strengths. In this case, a correction for H/l 
emission was determined from the measured [OIII] emission 
strength, as described earlier. We therefore apply the same 
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Figure B2. Comparison with the indices of Gonzalez (1993, 
G93). In both cases, H/3 is corrected for emission by adding 0.7 
times the [OIII] emission. The G93 values were measured from 
one or two long-slit position angles. The A3D values were mea¬ 
sured on circular apertures using the radii given in table 4.7 of 
G93, which give the luminosity-weighted circular aperture radius 
equivalent to the long-slit apertures. Errors, dashed lines and plot 
titles are as per Figure Bl. 


correction procedure, adding 0.7x [OIII] as measured using 
the GANDALF code to the measured H/3 index, with no 
additional emission correction. 

Instead of comparing Fe5015 (which is not in the Gon¬ 
zalez table), we compare Fe5270, which has the advantage 
of being free from emission corrections. To avoid the filter 
truncation effects described earlier, we convert our mea¬ 
sured Fe5270s values to the Lick Fe5270 index following 
Kuntschner et al. (2006): Fe5270 = 1.28xFe5270s -I- 0.03. 

The comparison of all three comparable indices is rather 
good, with the mean difference in all three indices consistent 
with zero within the scatter. The effective reduced val¬ 
ues are also close to unity, implying that the assumed errors 
are representative. The higher degree of consistency with 
Gonzalez (1993) over that with Trager et al. (1998) likely 
reflects more the more accurate aperture matching with the 
latter study, rather than some intrinsic difference in the mea¬ 
surement quality (with respect to the errors). 



iog,o(tJj 

Figure B3. Simulation of the effect of a fixed 3” aperture on the 
slope of the relation between metallicity and velocity dispersion. 
Green points show the i7e/2 values presented in Figure 8. The red 
points were measured assuming a fixed aperture, scaled to mimic 
the redshift of the MOSES sample Thomas et al. (2010). The effect 
of the fixed aperture is to include more of the lower metallicity 
regions in the outskirts of smaller, lower dispersion galaxies, thus 
steepening the slope compared to our values. Red and green lines 
show robust fits to the points of corresponding colour. The black 
dashed line is the relation from Thomas et al. (2010). The slopes 
match within the uncertainties, but the offsets differ by 0.2 dex. 


B2 Comparison with MOSES sample 

In Section 4.1 we compared our relations of SSP parameters 
and velocity dispersion to those of Thomas et al. (2010), and 
found that their metallicity relation was significantly steeper 
than ours. We can understand the difference in the metallic¬ 
ity trend by considering that no corrections were made by 
Thomas et al. (2010) for the fixed 3” diameter fiber used. 
In Figure B3 we simulate the effect of measuring metallic¬ 
ity for our sample with a 3” fiber at a redshift matching 
the mean redshift of the MOSES study (0.055). We overplot 
robust linear fits, excluding objects younger than 3.2 Gyr, 
which accounts for the redshift difference of our sample. The 
slopes agree within the errors, but a 0.2 dex offset in the re¬ 
lations remains. As in Kuntschner et al. (2006), we find that 
applying the models of either Schiavon (2007) or Thomas 
et al. (2005) yield very similar results, thus cannot explain 
the large metallicity offset (indeed, the offset becomes even 
slightly larger). Most likely it arises from the different com¬ 
bination of indices used by the two studies, but it is beyond 
the scope of this paper to explore this further. 


APPENDIX C: GENERIC APERTURE 
CORRECTIONS FROM THE ATLAS®° SAMPLE 

Since we have the privilege of extensive two-dimensional in¬ 
formation for a complete sample of galaxies, it is useful to 
investigate aperture corrections more generally, for use by 
studies without access to this spatial information. Figure Cl 
presents index growth curves as a function of Re for all ob¬ 
jects with maximal radial coverage of Re/2 or greater. The 
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Figure Cl. Curves of growth for each of our line indices, mea¬ 
sured within circular apertures of increasing radius, plotted as a 
function of the effective radius. The curves are normalized by the 
aperture value at i?e/2. Colours indicate the velocity dispersion 
within the normalization radius, ordered low (blue) to high (red). 
The central 1.2" are excluded, as are aperture sizes that are less 
than 85% filled with spectra. 


curves are normalized by the value at i?e/2. The colour of 
the curves are ordered in velocity dispersion from low (blue) 
to high (red). We do not include Fe5270 in this analysis due 
to the limited field coverage issue mentioned above. 

From this one can see differences between the red and 
blue curves, especially in the case of H/3 and Mgb. Com¬ 
paring this figure to the almost equivalent Figure 16 of 
Kuntschner et al. (2006), we see that there are many more 
galaxies with rising central H/3 strengths in our larger sam¬ 
ple, and these galaxies tend to have lower velocity disper¬ 
sions, although for a given velocity dispersion, there is a 
wide variety of slopes. 

We attribute this behaviour to the presence of a 
centrally-concentrated population of young stars that give 
rise to a steep gradient in H/3. If this is not taken into ac¬ 
count in the aperture correction, the strength of H/3 may be 
over-estimated. Figure C2 plots the curve-of-growth slopes 
(measured from robust line hts to the coloured curves in Fig¬ 
ure Cl) for all galaxies as a function of the H/3 line-strength 
measured within Re/8, for each of our three indices. There is 
clearly a relationship between the slopes of H/3 and Mg b with 
the Re/8 H/3 line-strength, presumably reflecting the known 



Figure C2. The power-law slope values for fits to all the galaxies 
in Figure Cl, plotted as a function of H/3 line-strength within an 
Re/S aperture. This illustrates how the size of aperture correction 
changes as a function of the central H/3 line-strength, due to the 
typically centrally-concentrated nature of the young stars. Dotted 
lines show the fixed corrections from Kuntschner et al. (2006). 


Table Cl. Coefficients of linear fits to aperture correction factors 
as a function of central H/3 strength, for use with equation C2. 


Index 

a 

b 

H/3 

-0.052 

0.069 

Fe5015 

0.021 

-0.105 

Mgb 

0.047 

-0.151 


age-sensitivity of these two lines. Thus, instead of providing 
a single aperture correction value for our three indices, we 
construct robust linear fits to the data in Figure C2, giving 
the corrections as a function of the H/3 lines-strength within 
an Re/8 aperture. Table Cl gives the coefficients of these lin¬ 
ear fits, which can be used to derive the aperture correction 
for each index such that: 

IndeXH/IndeX„orm = iR/RnormR (Cl) 

where: 

7 = aHPne/s + b (C2) 

and a and b are given in Table Cl. The H/3 dependence is 
important for central H/3 values larger than around 2 A, and 
neglecting the effect in this regime will lead to over-estimated 
Balmer line strengths, and under-estimated Mg 6. Below this 
threshold, however, the dependence is marginal, and adopt¬ 
ing a constant correction (such as that of Kuntschner et al. 
(2006), shown as dotted lines) gives similar results. 
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[MgFe50] 

Figure Dl. Large open diamonds show the Schiavon (2007) SSP 
model index predictions for a fixed age of 13.8 Gyr, as a func¬ 
tion of metallicity (increasing left to right) and abundance ratio 
(indicated by colour). The solid black line is a cubic fit to these 
predictions which describes the H/3 at this age as a function of 
[MgFeSOjL This line is used to find the offset in H/3 between the 
measured points (small diamonds) and the canonical ages of the 
Universe, as shown in Figure D2. For comparison, we also show the 
equivalent predictions from the models of Thomas et al. (2011), 
using the stellar evolutionary tracks of Cassis! et al. (1997, dot¬ 
ted line) and the Padova group (Girardi et al. 2000, dashed line). 
The mean measurement error is also indicated in the upper-right 
corner. 


APPENDIX D: GALAXIES OLDER THAN THE 
UNIVERSE? 

A number of galaxies show H/3 absorption values that imply 
ages older than the currently accepted age of the Universe 
(13.798 ± 0.037Gyr; Planck Collaboration, Ade et al. 2013). 
Here we ask the question: Are these objects consistent with 
this upper limit given our observational uncertainties? Even 
though the SSP models make predictions up to 17 Gyr, some 
data points still lie just outside the model grid. For this 
reason, we determine instead the differences of the measured 
H/3 to the H/3 predicted by the stellar population models for 
the age of the universe. 

We note that and age are not completely orthogo¬ 
nal, and there is a small but significant metallicity depen¬ 
dence. Abundance ratio dependences are small for the Schi¬ 
avon (2007) models used here. Figure Dl shows the line- 
strength predictions at the WMAP age for all metallicities 
and abundance ratios in the model grid (coloured diamonds). 
The solid line shows a cubic fit to the model points, giving 
the mean predicted H/3 value as a function of [MgFeSO]'. For 
comparison, we also show equivalent cubic fits derived from 
the “TMJ” models of Thomas et al. (2011) using the stellar 
evolutionary tracks of Cassisi et al. (1997, “TMJ-|-Cassisi”) 
and the Padova group (Girardi et al. 2000, “TMJ-|-Padova”). 
We shift all model predictions to the Lick system using the 
offsets from Schiavon (2007). This comparison shows that the 
data extend below predictions for the fiducial age for all mod¬ 
els, though the issue is somewhat less for the TMJ-|-Padova 
models, and worse for the TMJ-|-Cassisi models. 

Figure D2 shows the distribution of differences between 



Figure D2. Histogram showing the distribution of differences 
between the measured H/3 values and the model predictions for 
the canonical age of the Universe (13.8Gyr). The red line shows 
a Gaussian fitted to the lower portion of the distribution (fitted 
region is indicated), with a mean close to zero, and a dispersion 
similar to the observational error in H/3. From this we conclude 
that the galaxies appearing apparently older than the Universe 
are consistent with our measurement uncertainties. 


the observed H/3 values and the prediction for 13.8 Gyr from 
the Schiavon (2007) models used in this paper. Galaxies with 
negative differences have implied ages older than the Uni¬ 
verse. We then fit a Gaussian to the lower portion of this 
distribution, and find that the tail of negative points can 
be described by objects with a slightly positive mean off¬ 
set (i.e. younger than 13.8 Gyr), and a dispersion of 0.137A. 
This is quite consistent with our estimate of the H/3 measure¬ 
ment error, which is assumed to have a minimum of 0.1 A. 
Part of the additional observed dispersion may arise from 
uncertainties in the emission-line correction, in the form of 
emission lines below our detection limit, which is typically a 
few hundredths of an Angstrom at the signal-to-noise ratios 
considered here (Sarzi et al. 2006). We therefore conclude 
that the distribution of ages we measure are consistent with 
the fiducial age of the Universe after taking into account our 
measurement uncertainties, as well as systematic uncertain¬ 
ties in the models. 
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